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Abstract 

Q-( We present an analysis of the loop-induced magnetic dipole operator in the Randall- 

Q^ Sundrum model of a warped extra dimension with anarchic bulk fermions and an IR 

brane-localized Higgs. These operators are finite at one-loop order and we explicitly 
calculate the branching ratio for /i — )■ 67 using the mixed position/momentum space 
^ formahsm. The particular bound on the anarchic Yukawa and Kaluza-Klein (KK) 

scales can depend on the flavor structure of the anarchic matrices. It is possible for a 
^ generic model to either be ruled out or unaffected by these bounds without any fine 

O tuning. We quantify how these models realize this surprising behavior. We also review 

^ tree-level lepton fiavor bounds in these models and show that these are are in mild 

^ tension with the yU — )■ 67 bounds from typical models with a 3 TeV Kaluza-Klein scale. 

O Further, we illuminate the nature of the one-loop finiteness of these diagrams and 

^ show how to accurately determine the degree of divergence of a five-dimensional loop 

_ ^ diagram using both the five-dimensional and KK formalism. This power counting can 

^ be obfuscated in the four-dimensional Kaluza-Klein formalism and we explicitly point 

^ out subtleties that ensure that the two formalisms agree. Finally, we remark on the 

existence of a perturbative regime in which these one-loop results give the dominant 
contribution. 



Contents 



1 Introduction 

2 Review of anarchic Randall-Sundrum models 

3 Tree-level constraints from yU — )■ 3e and /i —t- e conversion 

3.1 Minimal RS model 

3.2 Custodially protected model 

4 Operator analysis of /i — )■ 67 

5 Calculation of /i — )■ 67 in a warped extra dimension 

5.1 Calculation of a 

5.2 Calculation of 6 

5.3 Modifications from custodial modes 

5.4 Constraints and tension 

6 Power counting and finiteness 

6.1 4D and 5D theories of bulk fields 

6.2 Bulk fields in the 5D formalism 

6.3 Bulk fields in the KK formalism 

6.4 Brane fields in the 5D formalism 

6.5 Brane fields in the KK formalism 

6.6 Matching KK and loop cutoffs 

6.7 Two-loop structure 

7 Outlook and Conclusion 

A Matching 5D amplitudes to 4D EFTs 

B Estimating the size of each diagram 

B.l Relative sizes of couplings 

B. 2 Suppression mechanisms in diagrams 

B. 3 Dimensionless integrals 

B. 4 Robustness against alignment 

C Analytic expressions 

Cl Dominant diagrams 

C. 2 Subdominant a coefficient diagrams 

C. 3 Subdominant h coefficient diagrams 

C. 4 Custodial Models 

D Position, momentum, and position/momentum space 
E Bulk Feynman Rules 

F Derivation of fermion propagators 

F.l Flat 5D fermion propagator 

F. 2 Warped 5D fermion propagator 

F. 3 Euclidean warped 5D fermion propagator 

G Finiteness of the brane-localized neutral Higgs diagram 

G. l Momentum space 

G. 2 Position/momentum space 



1 



1 Introduction 



The Randall-Sundrum (RS) set up for a warped extra dimension is a novel framework for models of 
electroweak symmetry breaking |1] . When fermion and gauge fields are allowed to propagate in the 
bulk, these models can also explain the fermion mass spectrum through the split fermion proposal 
[2]. In these anarchic flavor models each element of the Yukawa matrices can take natural 0{1) 
values because the hierarchy of the fermion masses is generated by the exponential localization of 
the fermion wave functions away from the Higgs field |3l H]. 

The same small wavefunction overlap that yields the fermion mass spectrum also gives hierar- 
chical mixing angles [21 E] and suppresses tree-level flavor-changing neutral currents (FCNCs) by 
the RS-GIM mechanism [21 E]- This built-in protection, however, may not always be sufficient to 
completely protect against the most dangerous types of experimental FCNC constraints. In the 
quark sector, for example, the exchange of Kaluza-Klein (KK) gluons induces left-right operators 
that contribute to CP violation in kaons and result in generic bounds of (9(10 — 20 TeV) for the 
KK gluon mass [Gj [TJ [H [9] . To reduce this bound one must either introduce additional structure 
(such as horizontal symmetries [TU] or flavor alignment [TT]) or alternately gain several 0{1) fac- 
tors [12] by promoting the Higgs to a bulk field, inducing loop-level QCD matching, etc. This 
latter approach is limited by tension with loop-induced flavor- violating effects [13|. 

The leptonic sector of the anarchic model is similarly bounded by FCNCs. Agashe, Blechman 
and Petriello recently studied the two dominant constraints in the lepton sector: the loop-induced 
/i — )■ 67 photon penguin from Higgs exchange and the tree-level contribution to /i — ?■ 3e and /i — )■ e 
conversion from the exchange of the Z boson KK tower [13] . These processes set complementary 
bounds due to their complementary dependence on the overall magnitude of the anarchic Yukawa 
coupling, Y^:. While /z — )■ 67 is proportional to due to two Yukawa couplings and a chirahty- 
fiipping mass insertion, the dominant contribution to /i — )■ 3e and /i — )■ e conversion comes from 
the nonuniversality of the Z boson near the IR brane. In order to maintain the observed mass 
spectrum, increasing the Yukawa coupling pushes the bulk fermion profiles away from the IR 
brane and hence away from the fiavor-changing part of the Z. This reduces the effective four- 
dimensional (4D) FCNC coupling so that these processes are proportional to Y~^. For a given KK 
gauge boson mass, these processes then set an upper and lower bound on the Yukawa coupling 
which are usually mutually exclusive. 

A key feature of the lepton sector is that one expects large mixing angles rather than the 
hierarchical angles in the Cabbibo-Kobayashi-Maskawa (CKM) matrix. One way to obtain this 
is by using a global flavor symmetry for the lepton sector [15] (see also [16]). Including these 
additional global symmetries can relax the tension between the two bounds. For example, imposing 
an A4 symmetry on the leptonic sector completely removes the tree-level constraints [15] . Another 
interesting possibility for obtaining large lepton mixing angles is to have the wavefunction overlap 
for the neutrino Yukawa peak near the UV brane [17] . For generic models with anarchic fermions, 
however, |14j found that the tension between /x — )■ 67 and tree-level processes (/i — )■ 3e and fi ^ e 
conversion) push the gauge boson KK scale to be on the order of 5-10 TeV. 

The main goal of this paper is to present a detailed one-loop calculation of the /i — )■ 67 penguin 
in the RS model with a brane-localized Higgs and to show that this amplitude is finite. 

To perform the calculation and obtain a numerical result we choose to work in the five- 
dimensional (5D) mixed position/momentum space formalism [181 UH]- This setup is natural for 
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calculating processes on an interval with brane-localized terms, as shown in Fig. [TJ In particular, 
there are no sums over KK modes, the chiral boundary conditions are fuUy incorporated in the 5D 
propagators, and the UV behavior is clear upon Wick rotation where the basis of Bessel functions 
becomes exponentials in the 4D loop momentum. The physical result is, of course, independent 
of whether the calculation was done in 5D or in 4D via a KK decomposition. We show exphcit 
one-loop finiteness in the KK decomposed theory and remark upon the importance of taking into 
account the correct number of KK modes relative to the momentum cutoff when calculating finite 
5D loops. 




Figure 1: A contribution to /i — )■ 67 from a brane-localized Higgs. The dashed line represents the 
Higgs while the cross represents a Yukawa coupling with a Higgs vev. 



The paper is organized as follows: We begin in Sections |2] and [3] by reviewing the flavor 
structure of anarchic Randall-Sundrum models and summarizing tree-level constraints on the 
anarchic Yukawa scale. We then proceed the analysis of /i — >■ 67. The dipole operators involved in 
this process are discussed in Section |4] and the relevant coefficient is calculated using 5D methods 
in Section [5j In Section [6] we discuss the origin of finiteness in these operators in both the 5D 
and 4D frameworks. We remark on subtleties in counting the superficial degree of divergence, the 
matching of the number of KK modes with any effective 4D momentum cutoff, and remark on 
the expected two-loop degree of divergence. We conclude with an outlook for further directions 
in Section [7} In Appendix |A] we highlight the matching of local 4D effective operators to nonlocal 
5D amplitudes. Next in Appendices |B] and [C] we give estimates for the size of each diagram and 
analytic expressions for the (next-to)leading /i — )■ 67 diagrams. Appendices D, E| and|F]focus on 
the formalism of quantum field theory in mixed position/momentum space, respectively focusing 
on a discussion of power counting, a summary of RS Feynman rules, and details on the derivation 
of the bulk fermion propagators. Finally, in Appendix [G] we explicitly demonstrate a subtle 
cancellation in the single-mass insertion neutral Higgs diagram that is referenced in Section [6] 



2 Review of anarchic Randall-Sundrum models 

We now summarize the main results for anarchic RS models. For a review see, e. g. Refs |20]. We 
consider a 5D warped interval z G [R, R'] with a U V brane aX z = R and an IR brane aX z = R! . 
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The metric is 



ds^ = (^^^ {dxi,dx,ri^"' - dz^), (2.1) 

where we see that R is also the AdS curvature scale so that R/R' ~ TcV/Mpi. Thesc conformal 
coordinates are natural in the context of the AdS/CFT correspondence but differ from the classical 
RS conventions z — Rex.p{ky) and k — l/R. The relevant scales have magnitudes R~^ ~ Mpi 
and R'~^ ~ TeV. Fermions are bulk Dirac fields which propagate in the fuU 5D space and can be 
decomposed into left- and right-handed Weyl spinors x V' via 

In order to obtain a chiral zero mode spectrum, these fields are subject to the chiral (orbifold) 
boundary conditions 

M^", R) = M^'. R') = o xr{x'', R) = Xi^(^^ R') = o, (2.3) 

where the subscripts L and R denote the SU{2)l doublet (L) and singlet (R) representations, i.e. 
the chirality of the zero mode. The fermion bulk masses are given hy c/R where c is a dimensionless 
parameter controlling the localization of the normalized 5D zero mode profiles, 

Xf (^,^) = ^ (D' {^y^fcX^c\x) and ^l^^f\x,z)=x^^l{x,z), (2.4) 

where we have defined the usual RS flavor function 



\l l-(R/R'y-2c- (2-5) 



We assume that the Higgs is locahzed on the IR brane. The Yukawa coupling is 



'Yuk 



= / (^) Ei{RYij)Lj- H + h.c. 



(2.6) 



where Yij is a dimensionless 3x3 matrix such that {Y^)ij = RYij is the dimensionful parameter 
appearing in the 5D Lagrangian. In the anarchic approach Y is assumcd to bc a random matrix 
with average elements of order Y^. After including all warp factors and rescaling to canonical 
fields the effective 4D Yukawa and mass matrices for the zero modes are 

so that the fermion mass hicrarchy is set by the /1^/2^ /3 structure for both left- and 
right-handed zero modes. In other words, the choice of c for each fermion family introduces 
additional fiavor structure into the theory which generates the zero mode spectrum while allowing 
the fundamental Yukawa parameters to be anarchic. 
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In the Standard Model the diagonalization of the fermion masses transmits the flavor structure 
of the Yukawa sector to the kinetic terms via the CKM matrix where it is manifested in the 
flavor-changing charged current through the boson. We shall use the analogous mass basis 
in Section [3] for our calculation of the Yukawa constraints from /i — )■ 3e and /x — )■ e conversion 
operators. The key point is that in the gauge basis the interaction of the neutral gauge bosons is 
flavor diagonal but not flavor universal. The different fermion wave functions cause the overlap 
integrals to depend on the bulk mass parameters. Once we rotate into the mass eigenbasis we 
obtain flavor changing couplings for the neutral KK gauge bosons. 

In the lepton sector this does not occur for the zero mode photon since its wavefunction remains 
flat after electroweak symmetry breaking and hence /i — )■ 67 remains a loop-level process. Thus for 
the primary analysis of this paper we choose a basis where the 5D fields are diagonal with respect 
to the bulk masses while the Yukawas are completely general. In this basis all of the relevant 
flavor-changing effects occur due to the Yukawa structure of the theory with no contributions from 
W loops. In the Standard Model, this corresponds to the basis before diagonalizing the fermion 
masses so that all flavor-changing effects occur through off-diagonal elements in the Yukawa matrix 
manifested as mass insertions or Higgs interactions. This basis is particularly helpful in the 5D 
mixed position/momentum space framework since the Higgs is attached to the IR brane, which 
simplifies loop integrals. 



3 Tree-level constraints from /i — > 3e and /i — )► e conversion 

For a fixed KK gauge boson mass Mkk, limits on /i — 3e and fi e conversion in nuclei 
provide the strongest lower bounds on the anarchic Yukawa scale K. These tree-level processes 
are parameterized by Fermi operators generated by Z and Z' exchange, where the prime indicates 
the KK mode in the mass basis. The effective Lagrangian for these lepton flavor-violating Fermi 
operators are traditionally parameterized as |2T] 



G 

+g^i{eLl^^lL){eB^^,eR)\ + -^e-i^{v - 075)// ^ gT^l^^ - ^''Ts)?, (3.1) 

where we have only introduced the terms that are non-vanishing in the RS set up, and use 
the normalization where = T3 — 2Q'^ sin^ 6. The axial coupling to quarks, a"^, vanishes in the 
dominant contribution coming from coherent scattering off the nucleus. The (73,4,5,6 are responsible 
for /i — )■ 3e decay, while the f, a are responsible for /i — )■ e conversion in nuclei. The rates are 
given by (with the conversion rate normalized to the muon capture rate): 

Br(/i ^ 3e) = 2{gl + gl) + gl + gl , (3.2) 

PeEeG%F^mla^Zfff „ „ „ 
Br(/. ^ e) = ' '-^QW + «'), (3.3) 

TT ^ i capt 

where the parameters for the conversion depend on the nucleus and are calculated in the Feinberg- 
Weinberg approximation [22] and we write the charge for a nucleus with atomic number Z and 
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neutron number as 



Qn = v''{2Z + N)+v\2N + Z). (3.4) 
. The most sensitive experimental constraint comes from muon conversion in '^Ti, for which 

10^ 

Ee Pe ~ m^, Fp ~ 0.55, Z^s ~ 17.61, Tcapt ~ 2.6 ■ — . (3.5) 

S 

We now consider these constraints for a minimal model (where /e^ = /e^, /^j, = Z^^.) and for a 
model with custodial protection. 

3.1 Minimal RS model 



In order to calculate the coefficients in the effective Lagrangian (3.1), we need to estimate the 
flavor-violating couplings of the neutral gauge bosons in the theory. In the basis of physical 
KK States all lepton flavor-violating couplings are the consequence of the non-uniformity of the 
gauge boson wave functions. Let us first consider the effect of the ordinary Z boson, whose wave 
function is approximately (we use the approximation (2.19) of |23] with a prefactor for canonical 
normalization) 



/i(°)(^) = ^ 



i? log f 



l + ^.^(l-21og-) 



(3.6) 



The coupling of the Z to fermions can be calculated by performing the overlap integral with the 



fermion profiles in (2.4) and is found to be 



After rotating the fields to the mass eigenbasis we find that the off-diagonal coupling of the Z 
boson to charged leptons is given by the nonuniversal term and is approximately 

9l,r ~ {dsM} ^'^eil = {9sm) 2(3 — 2c) ^ f'^L,Rffj.L,R- (3-8) 



Using these couplings one can estimate the coefficients of the 4-Fermi operators in (3.1) 



^3,4 = 2(7i,^A(°) (75,6 = 2^?l^7/?AW {v±a)=2gL,RA^^J, (3.9) 

where the gL,R are proportional to the left- and right-handed charged lepton couplings to the Z 
in the Standard Model, gi = —^ + s'^y and gn = s^. The Z' exchange contribution to /i — )■ 3e 
{fj, — )■ e) is a 15% (5%) correction and the 7' exchange diagram is an additional 5% (1%) correction; 
we shall ignore both here. We make the simplifying assumption that /e^ = /e^ and /^^ = /^^ 
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and then express these in terms of the Standard Model Yukawa couplings as / = ^J\/Y^. The 
expressions for the lepton flavor-violating processes are then 

Br(f, ^ e)T, = 2 . 10-1^ y fl-V . (3.11) 



The current experimental bounds are Br(/x — )• 3e) < 10 [21] and Br(/i — )■ e)Ti < 6.1 ■ 10 
|25j so that /i — e conversion provides the most stringent constraint, 



3 TeV\^ / 2 



< 0.5. (3.12) 



Mkk / VF* 

For a 3 TeV Z', the anarchic Yukawa scale must satisfy F* > 3.7, which agrees with [H] . 
3.2 Custodially protected model 



Since the bound in (3.12) is model dependent, one might consider weakening this constraint by 



having the leptons transform under the custodial group 

SU(2)^ X SU(2)^ X U(l)^ X Vlr, (3.13) 

where Pl/j is a discrete L -H- R exchange symmetry. Such a custodial protection was introduced 
in [26] to eliminate large corrections to the Zbb vertex in the quark sector. It was later found 
that this symmetry also eliminates some of the FCNCs in the Z sector [8] so that one might 
also expect it to alleviate the lepton flavor violation bounds. We shall now estimate the extent 
to which custodial symmetry can relax the bound on F^,. Further discussion including neutrino 
mixing can b e found in |27j . 

To custodially protect the charged leptons one choses the (L, R)x representation (2, 2)o for the 
left-handed leptons, (3, l)o© (1, 3)o for the charged right-handed leptons, and (1, l)o for the right- 
handed neutrinos. There are two neutral zero mode gauge bosons, the Standard Model Z and 7, 
and three neutral KK excitations, 7', Z' and Z^, where the latter two are linear combinations of 
the Z and Zx boson modes. The coupling of the left handed leptons to the ordinary Z and the 
Z' are protected since those couplings are exactly flavor universal in the limit where Plr is exact. 
The breaking of Plr on the UV brane leads to small residual contributions which we neglect. The 
remaining flavor-violating couplings for the left-handed leptons come from the exchange of Z^ 
and the 7', while the right-handed leptons are unprotected. 

Since {v — a) couples to right-handed leptons its coupling is unprotected and is the same as in 



(3.9). For (v + a), on the other hand, the leading-order effect comes from the Z^^^ component of 



the Zh, whose composition in terms of gauge KK states is [8] 

Zh = cos ^Z^^^ + sin ^Z^^^ + /3Z^°^ , (3.14) 



where is the flat zero mode Z-boson which does not contribute to FCNCs, cos ^ ~ y | ~ s^^/cw, 
and /3 is a small correction of order (9(f^/M^j^). The flavor-changing coupling of the KK gauge 
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bosons is analogous to that of KK gluons in [U], 



ati'' - (^Im)''''aJ;^(^) ^ fe)''Vlogf 7c/e,«/.,«, (3.15) 



V2 f' , i_2c w ^ . . ^ 0.7X1 



where 

7c = I dx x'-''^Jiixix) ^ ^^^^^^^^^ (3.16) 

Jiixi)Jo ^ ' Ji(xi)2(3-2c) ^ ' 

and xi = Mkk-R' is the first zero of Joix). The analogous 7*^^^ coupling is given by — )■ e. 
Taking into account the Zfj and 7''^'', the {v + a) effective couphng to left-handed leptons is 



(t; + a) = 2gL gKK ji^ \ cos 



e + ^cos^sinA A2^) + 24c^^,KK^^Af(^). (3.17) 



The cos^ sin^ term in the parenthesis represents the Z^'' component of the which couples to 
the quarks in the nucleus via 

Qfr = -^cw cose (5Z + 7N) - ^sw-{Z + N), gKK = (3.18) 

V 2 cos 4 g ^y\ogR'/R 

The gKK factor gives the universal (flavor-conserving) coupling of KK gauge bosons to zero mode 
fermions. is the electric charge of the nucleus normalized according to (3.3), = 2Z. 

Minimizing over the flavor factors /e^ ^ and /^^ ^ subject to the zero mode fermion mass 
spectrum and comparing to the experimental bound listed above (3.12 ), we find that the conversion 
rate must satisfy 



3 TeV\^ / 2 



Mkk j 

lowering the bound to K,, > 1 for a 3 TeV KK gauge boson scale. 



< 1.6. (3.19) 



4 Operator analysis of /i — > 67 

We work in 't Hooft-Feynman gauge (^ = 1) and a flavor basis where all bulk masses Cj are 
diagonal. The 5D amplitude for /i —t- 67 takes the form 

CH-L,a^'''E,FMN, (4.1) 

where it is understood that the 5D fields should be replaced by the appropriate external states 
which each carry an independent z position in the mixed position/momentum space formalism. 
These positions must be separately integrated over when matching to an effective 4D operator 



so that (4.1) can be thought of as a dimension-8 5D scattering amplitude whose prefactor C is 
a function of the external state positions, as explained in Appendix |A} When calculating this 
amplitude in the mixed position/momentum space formalism, the physical external state fields 
have definite KK number, which we take to be zero modes. The external field profiles and internal 
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propagators depend on 4D momenta and 2;-positions so that vertex z-positions are integrated from 
z = R to z = R' while loop momenta are integrated as usual. 

After plugging in the wave functions for the fermion and photon zero modes, including all 
warp factors, matching the gauge coupling, and expanding in Higgs-induced mass insertions, the 
leading order 4D operator and coefficients for — )■ 67 are 

^"yI^^/l. {aMY.kYi,Y,, + h,Y,j) I-eI^ a^"" Ef F^J + h.c. (4.2) 

The term proportional to three Yukawa matrices comes from the diagrams shown in Figs. |2] and 
|3} while the single- Yukawa term comes from those in Fig. |4} In the hmit where the bulk masses 
are universal, we may treat the Yukawas as spurions of the U (3)^ lepton flavor symmetry and note 
that these are the products of Yukawas required for a chirality-flipping, flavor- changing operator. 

In anarchic flavor models, however, the bulk masses for each fermion species is independent and 
introduce an additional flavor structure into the theory so that the U (3)^ lepton flavor symmetry 
is not restored even in the limit F — )■ 0. The indices on the dimensionless a^^ and hij coefficients 
encode this flavor structure as carried by the internal fermions of each diagram. Because the 
lepton hierarchy does not require very different bulk masses, both aku and h^j are nearly universal. 



Next note that the zero-mode mass matrix (2.7) introduces a preferred direction in flavor space 
which defines the mass basis. In fact, up to the non-universality of h^, the single- Yukawa term in 



(4.2) is proportional to — or aligned — with (2.7). Hence upon rotation to the mass basis, the off- 
diagonal elements of this term are typically much smaller than its value in the flavor basis [2H] and 
would be identically zero if the bulk masses were universal. Given a set of bulk mass parameters, 
the extent to which a specific off-diagonal element of the bij term is suppressed depends on the 
particular structure of the anarchic 5D Yukawa matrix. This is a novel feature since the structure 
of the underlying anarchic Yukawa is usually washed out in observables by the hierarchies in the 
/c flavor functions. 

On the other hand, a product of anarchic matrices typically indicates a very different direction 
in flavor space from the original matrix so that the Oj-,- term is not aligned and we may simplify 
the product to 

Y,(^kiY,uYlY,, = aV! (4.3) 

for each i and j. Here we have defined the prefactor a; different definitions can include an overall 
0{1) factor from the sum over anarchic matrix elements. We have used the anarchic limit and the 
assumption that neither ake nor bij vary greatly over realistic bulk mass values. This assumption 
is justified in Section |5] where we explicitly calculate these coefficients to leading order. Further, 
we have assumed that the scales of the anarchic electron and neutrino Yukawa matrices are the 
same so that (YE)ij ~ {YN)ij ~ K- 

To determine the physical fi ^ ej amplitude from this expression we must go to the Standard 
4D mass eigenbasis by performing a bi-unitary transformation to diagonalize the Standard Model 
Yukawa, 

= t/LA(^'"s)?7i, (4.4) 



9 



where the magnitudes of the elements of the unitary matrices Ul^fi are set, in the anarchic scenario, 
by the hierarchies in the flavor constants 

. /i. fe. < ,4.5, 

For future simpHcity, let us define the relevant part of the bijYij matrix after this rotation, 

= Y,iUL)2kbkMUl)n. (4.6) 

k,i 

The traditional parameterization for the /i — )■ 67 amphtude is written as [H] 

^^f^UL,Ra^''un,LF^,, (4.7) 



where ul,r are the left- and right-handed Dirac spinors for the leptons. Comparing (4.2) with 
(4.7) and using the magnitudes of the off-diagonal terms in the Ul rotation matrix in (4.5), we 
find that in the mass eigenbasis the coefficients are given by 

Cl = {aY^' + bV,) i?'2_|__^2m^^ J.^,, (4.8) 

Cn = {aY^' + 6K) i?'2_|__^2m^^ (4.9) 



The /i — )■ 67 branching fraction and its experimental bound are given by 

-1-11 



Br(/i ^ e7)thy = 77^^(|Cl|' + \Cr\'), (4.10) 

Br(;u^e7)exp< 1.2-10-". (4.11) 

While the generic expression for Br(/i — )■ 67) depends on the individual wave functions Jl-e, the 
product ClCr is fixed by the physical lepton masses and the relation C\ + C|. > 2ClCr so that 
one can put a lower bound on the branching ratio 

R'^V m, ,.0 ,,2/3 TeV^ ^ 



Br{f, -^e^)>6\aY,' + b ^—[ — ] ^ 5.1 ■ 10"« \aY,' + bn —— . (4.12) 
' ' Att \GfJ V ^KK 

Thus for a 3 TeV KK gauge boson scale we obtain an upper bound on F* 

3 TeV^ ^ 



Note that the b coefficient is independent of K,, so that sufficiently large b can rule out the assump- 
tion that the 5D Yukawa matrix can be completely anarchic — i. e. with no assumed underlying 
flavor structure — at a given KK scale no matter how small one picks Y^,. This is a new type of 
constraint on anarchic flavor models in a warped extra dimension. Conversely, if b is of the same 
order as a and has the opposite sign, then the bounds on the anarchic scale F* are alleviated. We 
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i/O, GO 



Figure 2: Neutral boson diagrams contributing to the a coefficient defined in (4.3). Fermion 



arrows denote the zero mode chirality, i. e. the SU(2) representation. External legs whose arrows 
do not point outward have an imphcit external mass insertion. Dotted hnes represent the fifth 
component of a bulk gauge field. Analytic forms for these diagrams are given in Appendix [C| 



will show below that b is typically suppressed relative to a but can, in principle, take a range of 
values between b = —0.5 and 0.5. For simphcity we may use the case 6 = O as a representative 
and plausible example, in which case the bound on the anarchic Yukawa scale is 



K < 0.12 lal 



(4.14) 



In Section |5.4| we quantify the extent to which the b term may affect this bound. Combined with 
the lower bounds on K,, from tree-level processes in Section [3| this bound typically introduces a 
tension in the preferred value of Y^, depending on the value of a. In other words, it can force one 
to either increase the KK scale or introduce additional symmetry structure into the 5D Yukawa 



matrices which can reduce a in (4.3) or force a cancellation in (4.13) 



5 Calculation of /i — 67 in a warped extra dimension 

In principle, there are a large number of diagrams contributing to the a and b coefficients even when 
only considering the leading terms in a mass insertion expansion. These are depicted in Figs. [2]-|4j 
Fortunately, many of these diagrams are naturally suppressed and the dominant contribution to 
each coefficient is given by the two diagrams shown in Fig. [5j Analytic expressions for the leading 
and next-to-leading diagrams are given in Appendix [C] along with an estimate of the size of each 
contribution. 

The flavor structure of the diagrams contributing to the b coefficient is aligned with the fermion 
zero-mode mass matrix [U [Hl |12]. The rotation of the external states to mass eigenstates thus 
suppresses these diagrams up to the bulk mass (c) dependence of internal propagators which point 
in a different direction in fiavor space and are not aligned. Since KK modes do not carry very 
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w'- 



Figure 3: Charged boson diagrams contributing to the a coefficient following the conventions in 
Fig. [2j Analytic forms for these diagrams are given in Appendix [C| 



strong bulk mass dependence, the diagrams which typically give the largest contribution after 
ahgnment are those which permit zero mode fermions in the loop. We provide a precise definition 
of the term "typically" in Section 5.2 



The Ward identity requires that the physical amplitude for a muon of momentum p to decay 
into a photon of polarization e and an electron of momentum p' takes the form 



M 



[{p + p'Y - {rrif, + nie)'^^] Up. 



(5.1) 



This is the combination of masses and momenta that gives the correct chirality-flipping tensor 



amplitude in (4.7). This simplifies the calculation of this process since one only has to identify 



the coefficient of the Up'{p + p')'^u term to determine the entire amplitude; all other terms are 
redundant by gauge invariance ^29j. The general strategy is to use the Chfford algebra and the 
equations of motion for the external spinors to determine this coefficient. This allows us to directly 
write the finite physical contribution to the amplitude without worrying about the regularization 



of potentially divergent terms which are not gauge invariant. In Section 6A we will further use 
this observation to explain the finiteness of this amplitude in 5D. 

In addition to the diagrams in Figs. |2}]4[ there are higher-order diagrams with an even number 
of additional mass insertions and brane-to-brane propagators. Following the Feynman rules in 
Appendix [E| each higher-order pair of mass insertions is suppressed by an additional factor of 



v 



0(10" 



(5.2) 



since we assume anarchic Yukawa matrices, K ~ 2. We are thus justified in considering only the 
leading-order terms in the mass insertion approximation. 

We now present the leading contributions to the a and b coefficients. Other diagrams give 
a correction on the order of 10% of these results. We provide explicit formulas and numerical 
estimates for the next-to-leading order corrections in Appendix O 
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Figure 4: Diagrams contributing to the b coefficient foUowing the conventions in Fig. |2| Not shown: 
zero mass-insertion diagram. Analytic forms for these diagrams are given in Appendix [Cj 



5.1 Calculation of a 

We now calculate the leading-order contribution to the amphtude to determine the a coefficient 
in (4.3). As discussed above, it is sufficient to compute the coefficient of the {p + p'Y term in 



the amphtude. The dominant contribution to a comes from the W boson diagrams in Fig. [5^. 
This is because diagrams with 5D gauge bosons are enhanced relative to the Higgs diagrams by 
a factor of InR'/R ~ 37. Further, the W diagrams are enhanced over the Z diagrams due to the 
size of their respective Standard Model couphngs to leptons. Additional suppression factors can 
arise from the structure of each diagram and are discussed in Appendix [B} Exphcit calculation 
confirms that the W loop with two internal mass insertions indeed gives the leading contribution 
to a. 

The charged and neutral boson diagrams have independent flavor structures, {YeY^Y]^)!^^ 
and (YeY^Ye) respectively. The anarchic Yukawa assumption implies that both of these terms 
should be of the same order, Y^. However one must remember that there may be a relative sign 
between these contributions depending on the specific anarchic y^r and Ye matrices. In other 
words, a = cichargcd ± ^neutral where the sign cannot be specified generically. However, because 
«neutral <^ «charged, we iguore the ueutral boson loops, though these neutral boson diagrams may 
become appreciable if one allows a hierarchy between the overall scales of the Yat and Ye matrices. 

The W loop in Fig. [5^ contains an implicit mass insertion on the external muon leg. As 
explained in Appendix [B| the 5D fermion propagator between this mass insertion and the loop 
vertex is dominated by the KK mode which changes fermion chirality. This is because the chirality- 
preserving piece of the propagator goes like p. Invoking the muon equation of motion gives a factor 

of fil'\vR')f^^ ~ {m^j_R') for the external leg. This is much smaller than the f'^\vR')f'l^^^ factor 
from the chirality-flipping part of the propagator. Compared to the mass insertion connecting 
the zero mode external muon to a KK intermediate state, the mass insertion connecting two zero 
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fa) W 



(b) W^,W 



Figure 5: The leading diagrams contributing to the a and h coefficients. following the same 
conventions as Fig. [2j 

mode fermions is smaller by a factor of the exponentially suppressed zero mode profil^ 
Using the Feynman rules in Appendix [El the amplitude this diagram is 



ip+p') iQ^2 J^L^^* J-CE, y/2\2 R 



71 



'2MW "p 



(5.3) 



where huiw = —0.31 is a dimensionless loop integral. Taking R'v/\/2 = .17 and \n{R' /R) = 

a = -0.065. (5.4) 



7.3, the a coefficient in (4.3) is 



5.2 Calculation of b 

As discussed above, the diagrams contributing to b are sensitive to the structure of the anarchic 
Yukawa matrix relative to that of the non-universal internal bulk fermion masses. For example, 
if the bulk mass parameters were universal, then the b coefficient operator would be aligned and 
the off-diagonal element would vanish. The sign of this off-diagonal term is a function of the 
initial anarchic matrix so that the b term may interfere constructively or destructively with the a 
term calculated above. We numerically generate anarchic matrices whose elements have random 
sign and random values between 0.5 and 2 to determine the distribution of probable Yukawa 
structures. Such a distribution is peaked about zero so that the choice 6 = O is a reasonable 
simplifying assumption. For a more detailed description of the range of bounds accessible by the 
anarchic RS scenario, one may use the la value of |6| as characteristic measure of how large an 
effect one should expect from generic anarchic Yukawas. 

The dominant contribution to the b coefficient has a and W propagating in the loop and 



is shown in Fig. [5p. Following the analysis in in Appendix |B.4[ this diagram can have zero mode 
fermions propagating in the loop and hence is particularly sensitive to the bulk mass parameters 
of the internal fermions being summed in the loop. This, in turn, implies that the diagram is 
more robust against alignment and can give an appreciable contribution after rotating to the zero 



^We thank Martin Beneke, Paramita Dey, and Jiirgen Rohrwild for pointing this out. 
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mode mass basis. The amplitude associated with this diagram is 

= ) {^J^^-^^^^^'^r^iP + PTu,, (5.5) 

where the fermion bulk mass dependence of the dimensionless integral is sketched in Fig. [6j The 
shape of this plot can be understood by recalling that the propagator has Dirichlet boundary 
conditions so that the loop has most of its support off of the IR brane. We shall approximate the 
fermion contribution by the zero mode only since this is the mode with the most dependence on 
c and hence the most resilient to alignment. The internal fermion gives its largest contribution 
when its profile has the biggest overlap with the rest of the loop integral. For c < 0.5 the zero 
mode profile leans towards the IR brane and overlaps with the rest of the loop. As c is taken above 
c ~ 0.5, the zero mode profile leans toward the UV brane and the overlap with the rest of the loop 
is exponentially suppressed with c; this explains the sharp drop in Fig. [6j In the minimal model, 
Ce lies on the ridge and so the integral gives a very different value for c = Cg versus c = c^^r- This 
difference gives a non-universal contribution to the b coefficient that is not aligned when passing 
to the zero mode mass basis. 



I(c) 




Figure 6: Plot of the dimensionless integral / for the W diagram that gives the leading 
contribution to the b coefficient after alignment. Black dots correspond to values for minimal 
model bulk masses. Because the electron has a large bulk mass Ce ~ 0.66, it falls the ridge and 
the value of the integral is different from /(c^) or I{cr), giving a large contribution an appreciable 
contribution after alignment. 

After scanning over anarchic matrices as defined above, the la value of the dimensionless 
integral is /g^ivi/s =0.2. This gives a Icr value for the b coefficient, 

jfei"! = 1.5. (5.6) 

Note that the large contribution of the W diagram to the la value of the b coefficient depended 
on the minimal model Cg value falling in a range where J(c) is sharply decreasing so that there is 
a large non-universal piece. The W diagram and one-mass-insertion Z diagrams also permit 
zero mode fermions in the loop and can also in principle give appreciable contributions but are 
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much smaller than the W diagram in the minimal model. In general, one may construct non- 
minimal models which allow cl^ to take values closer to the other left-handed bulk masses. In 
this way one could reduce the contribution of the W diagram increase the importance of the 
Z diagram. For our purposes we shall take the la value of the minimal model b coefficient as a 
representative benchmark for a plausible general estimate for the generically allowed range of b. 



5.3 Modifications from custodial modes 



In Section 3^ it was shown that custodial symmetry weakens the bounds from tree-level FCNCs. 
Since we would like to assess the tension between tree- and loop-level bounds, we should also 
examine the effect of the additional custodial modes on /i — )■ 67. These additional diagrams are 
described by the same topologies as those in Figs. |2]-[4] but differ by replacing internal lines with 
custodial bosons and fermions. The expression for the amplitude differs by coupling constants 
and the use of propagators with different boundary conditions, but not in the overall structure of 
each amplitude and so are straightforward to extract from the minimal model expressions. The 
leading topologies are unchanged so that it is sufficient to consider the custodial versions of the 
diagrams in Fig. [5j 

For the two-mass-insertion W diagram, there are two additional diagrams with custodial 
fermions: one with a Wl and the other with a Wr in the loop. The Plr symmetry enforces 
that the couplings are identical while the different boundary conditions modify the definitions of 
the internal propagators so that the only difference comes from the value of the dimensionless 



integral in (5.3). The each diagram contributes a dimensionless integral / = —0.2, so that the a 



coefficient is modified to 

Ocust. = -0.15. (5.7) 

For the b coefficient diagram we note that leptons carry no U(l)x charge so that any vertex 
between an external state fermion and a custodial boson must also have a custodial fermion. 
However, since custodial fermions are chosen to have no zero modes, the contribution of these 
diagrams is aligned with the light fermion mass matrix so that their contribution to b is negligible. 



5.4 Constraints and tension 



We can now estimate the upper bound on the anarchic Yukawa scale in (4.13), 

2 



\aY^ +b 



3 TeV 



Ml 



KK 



< 0.015. 



(4.13) 



First let us consider the scenario where the b coefficient takes its statistical mean value, 6 = 0, 
and Mkk = 3 TeV. In this case the minimal model suffers a (9(10) tension between the tree-level 
lower bound on F* and the loop-level upper bound. 



1; > 4 1; < 0.5. 

The custodial model slightly alleviates this tension, 

F* > 1.25 F < 0.3. 



(5.8) 
(5.9) 
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(a) Minimal model 



(b) Custodial model 



Figure 7: Bounds on the anarchic Yukawa and KK scales in the minimal (a) and custodial (b) 



models from tree- and loop-level constraints, (3.12), 



(3.19), and (4.13). Each curve rules out the 
The thick solid straight 



region to its left. The solid hyperbola is the appropriate tree- level bound. 
line is the 6 = loop-level bound. The red dashed (blue dotted) curve is the loop-level bounds in 
the case where b has the same (opposite) sign as a and takes its la magnitude |6| = \b\ia = 1.5. 



These discrepancies should be interpreted as an assessment on the extent to which the 5D Yukawa 
matrices may be generically anarchic. The tension in the bounds above imply that for Mkk = 
3 TeV, one must accept some mild tuning in the relative sizes of the 5D Yukawa matrix. This is 
shown by the hyperbola and solid line in Fig. [7} 

Alternately, one may ask that assuming totally anarchic Yukawas, what is the minimum value 
of Mkk for which the tension is alleviated? In the minimal model the tree- and loop-level bounds 
allow mutually consistent Yukawa scales for Mkk > 6 starting at Y = 1. Similarly, for the 
custodial model the tree- and loop-level bounds allow consistent values for Mkk > 4.75 starting 
at y = 0.5. 

Next one may consider the effect of the b coefficient which is sensitive to the particular flavor 
structure of the anarchic 5D Yukawa matrix relative to the choice of fermion bulk mass parameters. 
The la range of b values for randomly generated anarchic matrices is 6 G (—1.5, 1.5). Note that 
this can dwarf the a coefficient and the bound in (4.13). It is thus possible for statistically likely 
values of b of either sign to violate the bound independent of the value of Y^. Such a value of b 
would directly constrain the KK scale. For the la value this sets Mkk ^ 30 TeV; this can be 
seen from the intersection of the red dashed lines and blue dotted lines with the horizontal axes 
in Fig. [7) 

The most interesting range for b, however, is the statistically likely region near 6^0. In 
particular, it is possible that the flavor structure contributing to the b coefficient is such that the 
aligned diagram in Fig. [5]is of the same order of magnitude and opposite sign as the other diagram. 



17 



In this case one can relax the constraints on and Mkk- Of course it is equally likely that such 
a term would exacerbate these constraints so that for any specific choice of Yukawa matrices it 
is necessary to check the b term. For broad model-building purposes, however, the key point is 
that the flavor structure for each diagram in Fig. [5] is effectively independent and so one can find 
appreciable regions of flavor-space (spanned by the spurions Ye, Y^, and the bulk masses) where 
one avoids significant tension between the tree- and loop-level constraints. 

In Fig. [7| the red dashed hne shows the bound when b takes its Icr magnitude and has an 
opposite sign from a; the cusp at Mkk = O represents the case where the a and b terms canceL 
The blue dotted hne shows the case where b takes its la magnitude and has the same sign as 
a. What is important to note is that as one takes |6| less than these hnes continuously 

converge upon the straight hne corresponding to 6 = O so that any combination of F* and Mkk 
between the upper red dashed hne and the blue dotted hne can be plausibly achieved within 
the anarchic paradigm. Let us make the caveat that the above values are estimates at (9(10%) 
accuracy. Specific results depend on model-dependent factors such as the extent to which the 
matrices are anarchic, the relative scale of the charged lepton and neutrino anarchic values, or 
extreme values for bulk masses. For completeness we provide analytic formulas for the leading 
and next-to-leading order diagrams in Appendix [C} 

6 Power counting and finiteness 

We now develop an intuitive understanding of the finiteness of this 5D process, highlight some 
subtleties associated with the KK versus 5D calculation of the loop diagram^ and estimate the 
degree of divergence of the two-loop result. Our primary tool is nai've dimensional analysis, from 
which we may determine the superficial degree of divergence for a given 5D diagram. Special care 
is given to the treatment of brane-localized fields and the translation between the manifestly 5D 
and KK descriptions. 

6.1 4D and 5D theories of bulk fields 

It is instructive to review key properties of /i — )■ 67 in the Standard Model. This amplitude was 
calculated by several authors [221 El] • Two key features are relevant for finiteness: 

1. Gauge invariance cancels the leading order divergences. The Ward identity requires 
qfj,M.^ = O, where Ai^ is the amplitude with the photon polarization peeled off and 

is the photon momentum. This imposes a nontrivial g-dependence on M and reduces the 
superficial degree of divergence by one. 

2. Lorentz invariance prohibits divergences which are odd in the loop momentum, k. In other 
words, j d'^k^/k'^"' = 0. After accounting for the Ward identity, the leading contribution to 
the dipole operator is odd in k and thus must vanish. Specifically, one of the $ terms in a 
fermion propagator must be replaced by the fermion mass m. 

^The finiteness of dipole operators has been investigated in gauge-higgs unified models where a higher- 
dimensional gauge invariance can render these terms finite 30]. Here we do not assume the presence of such 
additional symmetries. 
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Recall that the chiral structure of this magnetic operator requires an explicit internal mass 
insertion. In the Standard Model this is related to both gauge and Lorentz invariance so that it 
does not give an additional reduction in the superficial degree of divergence. Before accounting for 
these two features, naiVe power counting in the loop integrals appears to suggest that the Standard 
Model amplitude is logarithmically divergent from diagrams with two internal fermions and a single 
internal boson. Instead, one finds that these protection mechanisms force the amplitude to go as 
where M is the characteristic loop momentum scale. 
We can now extrapolate to the case of a 5D theory. First suppose that the theory is modified to 
include a noncompact fifth dimension: then we could trivially carry our results from 4D momen- 
tum space to 5D except that there is an additional loop integral. By the previous analysis, this 
would give us an amplitude that goes as M^^ and is thus finite. Such a theory is not phenomeno- 
logically feasible but accurately reproduces the UV behavior of a bulk process in a compact extra 
dimension so long as we consider the UV limit where the loop momentum is much larger than 
the compactification and curvature scales. This is because the UV limit of the loop probes very 
small length scales that are insensitive to the compactification and any warping. This confirms 
the observation that /i — 67 in Randall-Sundrum models with all fields (including the Higgs) in 
the bulk is UV-finite [H] . In the case where there are brane-localized fields, this heuristic picture 
is complicated since the /i — ?• 67 loop is intrinsically localized near the brane and is sensitive to 
its physics; we address this issue below. 



6.2 Bulk fields in the 5D formalism 

We may formalize this power counting in the mixed position/momentum space formalism. This 
also generalizes the above argument to theories on a compact interval. Each loop carries an 
integral d^k and so contributes +4 to the superficial degree of divergence. We can now consider 
how various features of particular diagrams can render this finite. 



1. Gauge invariance {p + p')- As argued above and shown explicitly in (5.1), the Ward 
identity identifies the gauge invariant contribution to this process to be proportional to 
{p + p')'^, which reduces the overall degree of divergence by one. 

2. Bulk Propagators. The bulk fermion propagators in the mixed position/momentum space 
formalism have a momentum dependence of the form ^/k ^ 1 while the bulk boson propa- 
gators go like l/k. This matches the power counting from summing a tower of KKmodes. 



Note that this depends on k = y/k'^ so that the Lorentz invariance in Section QA_ for a 
noncompact extra dimension is no longer valid. 

Bulk vertices (dz), overall z-momentum conservation. Each bulk vertex carries an 
integral over the vertex position which brings down an inverse power of the momentum 
flowing through it. This can be seen from the form of the bulk propagators, which depend 
on z in the dimensionless combination kz up to overall warp factors. In the Wick-rotated UV 
limit, the integrands reduce to exponentials so that their integrals go like l/k. In momentum 
space this suppression is manifested as the momentum-conserving 6 function in the far UV 
limit where the loop momentum is much greater than the curvature scale. 
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An alternate and practical way to see the l/k scaling of an individual dz integral comes 
from the Jacobian as one shifts to dimensionless integration variables, 



y = kER' X = kEZ (6.1) 

so that y e [O, 00] plays the role of the loop integrand and x G [yR/R',y] plays the role of 
the integral over the interval extra dimension. These are the natural objects that appear as 
arguments in the Bessel functions contained in the bulk field propagators, as demonstrated 



in Appendix F. 3 In these variables each dx brings down a factor of l/y from the Jacobian 
of the integration measure. These variables are natural choices because they relate distance 
intervals in the extra dimension to the scales that are being probed by the loop process. The 
physically relevant distance scales are precisely these ratios. 

4. Overall 5;-momentum conservation. We must make one correction to the bulk vertex 
suppression due to overall 2;-momentum conservation. This is most easily seen in momen- 
tum space where one (5-function from the bulk vertices conserves overall external momentum 
in the extra dimension and hence does not affect the loop momentum. In mixed posi- 
tion/momentum space this is manifested as one dz integral bringing down an inverse power 
of only external momenta without any dependence on the loop momentum. We review this 
in Appendix|D| where we discuss the passage between position and momentum space. The 
overall 2;-momentum conserving 5-function thus adds one unit to the superficial degree of 
divergence to account for the previous overcounting of ~ l/Zc suppressions. 

5. Derivative coupling. The photon couples to charged bosons through a derivative coupling 
which is proportional to the momentum flowing through the vertex. This gives a contribution 
that is linear in the loop momentum, fc^. 

6. Chirality: mass insertion, equation of motion. To obtain the correct chiral structure 
for a dipole operator, each diagram must either have an explicit fermion mass insertion 
or must make use of the external fermion equation of motion (EOM). For a bulk Higgs 
field, each fermion mass insertion carries a dz integral which goes like l/k. As described 
in Section |5} the use of the EOM corresponds to an explicit external mass insertion. Thus 
fermion chirality reduces the degree of divergence by one unit. 

We may now straightforwardly count the powers of the loop momentum to determine the 
superficial degree of divergence for the case where the photon is emitted from a fermion (one 
boson and two fermions in the loop) or a boson (two bosons and one fermion in the loop). The 
latter case differs from the former in the number of boson propagators and the factor of k'^ in the 
photon Feynman rule. 
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Neutral Charged 





Boson 


Boson 


Loop integral (d'^k) 


+4 


+4 


Gauge invariance (p + p') 


-1 


-1 


Bulk fermion propagators 








Bulk boson propagator 


-1 


-2 


Bulk vertices (dz) 


-3 


-3 


Overall 2;-momentum 


+ 1 


+1 


Derivative coupling 





+1 


Mass insertion/EOM 


-1 


-1 


Total degree of divergence 


-1 


-1 



The WH"^ diagram in Fig. [4] is a special case since it has neither a derivative coupling nor an 
additional chirality flip, but these combine to make no net change to the superficial degree of 
divergence. We confirm our counting in Section 6.1 that the superficial degree of divergence for 
universal extra dimension where all fields propagate in the bulk is — 1 so that the fiavor-changing 
penguin is manifestly finite. 

Before moving on to the case of a brane-localized boson, let us remark that this bulk counting 
may straightforwardly be generalized to the case of a bulk boson with brane-localized mass inser- 
tions. To do this, we note that the brane-localized mass insertion breaks momentum conservation 
in the z direction and this no longer contributes -|-1 to the degree of divergence. On the other 
hand, each mass insertion no longer contributes —1 from the dz integral so that the changes in the 
"overall z-momentum" and "mass insertion/EOM" counting cancel out. We find that diagrams 
with a bulk gauge boson and brane-localized mass insertions have the same superficial degree of 
divergence as the lowest order diagrams in a bulk mass insertion expansion. 



6.3 Bulk fields in the KK formalism 

All of the power counting from the 5D position/momentum space formalism carries over directly 
to the KK formalism with powers of mKK treated as powers of k. The position/momentum space 
propagators already carry the Information about the entire KK tower as well as the profiles of 
each KK mode. Explicitly converting from a 5D propagator to a KK reduction, 

AMk,z,z') = Y,f-\z)A^^l{k)f-\z'), (6.2) 

n 

where /'•"^ is the profile of the rfi^ KK mode. The sum over KK modes is already accounted for in 
the 5D propagator; for example, for a boson A|^j^ ~ l/fc^ while ~ 1/fc. The vertices between 
KK modes are given by the dz integral over each profile, which reproduces the same counting 
since each profile depends on z as a function of m^i^z. Conservation of 2;-momentum is replaced 
by conservation of KK number in the UV limit of large KK number. 

Indeed, it is almost tautological that the KK and position/momentum space formalisms should 
match for bulk fields since the process of KK reducing a 5D theory implicitly passes through 
the position/momentum space construction. This will become slightly more nontrivial in the 
case of brane-localized fields. We shall postpone a discussion of mixing between KK states until 
Section 16.51 
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(a) Mass insertion before photon (b) Mass insertion after photon 



Figure 8: One-mass-insertion neutral scalar diagrams. The leading order /c-dependence of each 
diagram cancels when the two are summed together. 



6.4 Brane fields in the 5D formalism 

The power counting above appears to fail for loops containing a brane-locahzed Higgs field. The 
brane-locahzed Higgs propagator goes hke l/k"^ rather than Ijk for the bulk propagator, but this 
comes at the cost of two vertices that must also be brane-locahzed, thus negating the suppression 
from the dz integrals. The charged Higgs has two brane-localized Higgs propagators, but loses 
a third dz integral from the brane-localized photon emission. Finally, there are no additional 
contributions from the brane-localized fermion mass insertions nor are there any corrections from 
the conservation of overall z-momentum since it is manifestly violated by the brane-localized 
vertices (see Appendix |D] for a detailed discussion). In the absence of any additional brane effects, 
both types of loops would be logarithmically divergent, as discussed in p^. 

Fortunately, two such brane effects appear. First consider the two neutral Higgs diagrams in 
Fig. |2} The diagram with no mass insertion requires the use of an external fermion equation of 
motion which still reduces the superficial degree of divergence by one so that it is finite. The 
diagram with a single mass insertion is finite in the Standard Model due to a cancellation between 
the Higgs and neutral Goldstone diagrams, as discussed in Section [5j More generally, even for a 
single type of brane-localized field, there is a cancellation between diagrams in Fig. [8] where the 
photon is emitted before and after the mass insertion. This can be seen by writing down the Dirac 
structure coming from the fermion propagators to leading order in the loop momentum, 

Ma ~ ^Y^^-kYk^ = k^ih^-Yft) (6.3) 
Mb ~ ^^Y^-^kYk = k^Y^-h^) (6.4) 

The terms with three factors of ^ are contributions where "correct-chirality" fermions propagate 
into the bulk, while the terms with only one |^ are contributions where "wrong-chirality" fermions 
propagate into the bulk. The structure of the latter terms comes from the 'y^dz term in the 
Dirac operator. The structures above multiply scalar functions which, to leading order in k, are 



identical for each term. From the Clifford algebra it is clear that (6.3) and (6.4) cancel so that the 



contribution that is nonvanishing in the UV must be next-to-leading order in the loop momentum. 



In Appendix G this cancellation is connected to the chiral boundary conditions on the brane and is 
demonstrated with explicit flat-space fermion propagators. We thus find that the brane-localized 
neutral Higgs diagrams have an additional —1 contribution to the superficial degree of divergence. 
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Next we consider the charged Goldstone diagrams. These diagrams have an additional mo- 
mentum suppression coming from a positive power of the charged Goldstone mass appearing 
in the numerator due to a cancellation within each diagram. In fact, we have aheady seen in 



Section 5.1 how such a cancellation appears. For the single-mass-insertion charged Goldstone dia- 



gram in Fig. |3| we saw in (B.2) that the form of the 4D scalar propagators and the photon-scalar 
vertex cancels the leading-order loop momentum term multiplying the required {p + p')^. The 
cancellation introduces an additional factor of M^/i^k"^ — M^) so that the superficial degree of 
divergence is reduced by two. Note that the position/momentum space propagators for a bulk 
Higgs have a different form than that of the 4D brane-localized Higgs and do not display the same 



cancellation. In the KK picture this is the observation that the cancellation in (B.2) takes the 
form M^y^/lk"^ — M^j^), which does not provide any suppression for heavy KK Higgs modes. 

Finally, the diagrams where the photon emission vertex mixes the W and brane-localized 
charged Goldstone are special cases. The photon vertex carries neither a dz integral nor a k'^ 
Feynman rule and hence makes no net contribution to the degree of divergence. A straightforward 
counting including the brane-localized Goldstone, bulk W, and the single bulk vertex thus gives 
a degree of divergence of — 1. 

We summarize the power counting for a brane-localized Higgs as foUows: 
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Loop integral {d'^k) 


+4 
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-1 





-1 
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Total degree of divergence 


-1 
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-1 



It may seem odd that the brane-localized charged Higgs loop has a different superficial degree of 
divergence than the other 5D cases, which heretofore have all been —1. This, however, should not 
be surprising since the case of a brane-localized Higgs is manifestly different from the universal 
extra dimension scenario. It is useful to think of the brane-localized Higgs as a limiting form of 
a KK reduction where the zero mode profile is sharply peaked on the IR brane. The difference 
between the bulk and brane-localized scenarios corresponds to whether or not one includes the 
rest of the KK tower. 



6.5 Brane fields in the KK formalism 

Let us now see how the above power counting for the brane-localized Higgs manifests itself in the 
Kaluza-Klein picture [H]. Observe that this power counting for both the W-H^ and the charged 
boson loops are trivially identical to the 5D case due to the arguments in Section 6^ For example, 
the cancellation is independent of how one treats the bulk fields. The neutral Higgs loop. 
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however, is somewhat subtle since the "chiral cancellation" is not immediately obvious in the KK 
picture. 

We work in the mass basis where the fermion hne only carries a single KK sum (not independent 
sums for each mass insertion) and the zero mode photon couphng preserves KK number due to 
the flat A^^^ profile. In this basis the internal fermion line carries one KK sum and it is suflicient 
to show that for a single arbitrarily large KK mode the process scales like The four- 



dimensional power counting in Section 6.1 appears to give precisely this, except that Lorentz 
invariance no longer removes a degree of divergence. This is because this suppression came from 
the replacement of a loop momentum ^ by the fermion mass m. For an arbitrarily large KK 
mode, the fermion mass itself is the loop momentum scale and so does not reduce the degree of 
divergence. In the absence of any additional suppression coming from the mixing of KK modes, it 
would appear that the KK power counting only goes like 1/Mkk so that the sum over KK modes 
should be logarithmically divergent, in contradiction with the power counting for the same process 
in the 5D formalism. 

We shall now show that the pair of Yukawa couplings for the neutral Higgs also carries the 
expected l/k factor that renders these diagrams finite and allows the superficial degrees of diver- 
gence to match between the KK and 5D counting. It is instructive to begin by defining a basis 
for the zero and first KK modes in the weak (chiral) basis. We denote left (right) chiral fields of 
KK number a by Xlr (V'!"]?) where the L, R refers to SU(2)i doublets and singlets respectively. 
We can arrange these into vectors 

X = (xi°>. xs. xS) i- = *i?) . (6.5) 

where i runs over flavors. It is helpful to introduce a single index J = 3a + i where i = 1, 2, 3 
according to fiavor and a = 0, 1, 2 according to KK mode (writing a = 2 to mean the first KK 
mode with opposite chirality as the zero mode). Thus the external muon and electron are X2 and 
ipi respectively, while an internal KK mode takes the form xj ot ipj with J > 3. This convention 



in (6.5) differs from that typically used in the literature (e. g. P^) in the order of the last two 
elements of ip. This basis is useful because the KK terms are already diagonal in the mass matrix 
(^Mx + h.c.), 



(6.6) 










1 

m 


M = 1 




Mkk,i 











Mkk,2 



where each element is a 3 x 3 block in fiavor space and we have written 



m = -^Ayji^^ « Mkk, (6.7) 

with indices as appropriate and Mkk diagonal. Let us define e = v /Mkk to parameterize the 
hierarchies in the mass matrix. For a bulk Higgs, these terms are replaced by overlap integrals 
and the M32 block is nonzero, though this does not affect our argument. Note that Mkk.i and 
Mkk, 2 are typically not degenerate due to C (m) differences in the doublet and singlet bulk masses. 
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X2 J J ( ^ J J f. Vl 

Figure 9: The fermion line in the mass basis for diagrams with an internal KK mode (J > 3). For 
simplicity we do not show the internal photon insertion. 



In the gauge eigenbasis the Yukawa matrix is given by 



-M 
v 



















i) 


A/kk=0 


^0 








(6. 



where we have assumed fL,fR,y* ~ for simplicity since the hierarchies in the /*-°''s do not 
afFect our argument. The 1 elements thus refer to blocks of the same order of magnitude that are 
not generically diagonal. The O blocks must vanish by gauge invariance and chirality. 



We now rotate the fields in (6.5) to diagonalize the mass matrix (6.6); we indicate this by a 
caret, e. g. x- this basis the Yukawa matrix is also rotated y ^ y. The fermion line for this 
process is shown in Fig. [9j the Yukawa dependence of the amplitude is 

M ~ yijyj2- (6.9) 



First let us note that in the unrealistic case where y = y, one of the Yukawa factors in (6.9) is 
identically zero for all internal KK modes, J > 3. One might then expect that the mass rotation 
would induce a mixing of the zero modes with the KK modes that induces O {e) blocks into the 
Yukawa matrix, 



1 



(6.10) 



If this were the case then the product yijyj2 would not vanish, but would be proportional to 
e ~ 1/Mkkj which is precisely the KK dependence that we wanted to show. While this intuition 
is correct and captures the correct physics, the actual Yukawa matrix in the mass basis has the 
structure (c. f. (67) in [H]) 




1 + e -1 



(6.11) 



The new 0{1) elements come from the large rotations induced by the m^^ and m^^ blocks. These 
factors cancel out so that we still have the desired yijyj2 ~ e relation. Physically this is because 
these 0{1) factors come from the "large" rotation from chiral zero modes to light Dirac SM 
fermions. Thus they represent the "wrong-chirality" coupling of the external states induced by 
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the usual mixing of Weyl states from a Dirac mass. This does not include the mixing with the 
heavy KK modes, which indeed carries the above e factors so that the final result is 



yijyj2 ~ e ~ T-^, (6.12) 

giving the correct —1 contribution to the superficial degree of divergence for the neutral Higgs 
diagrams to render them manifestly finite. 

A few remarks are in order. First let us emphasize again that promoting the Higgs to a bulk 
field makes the 3-2 block of the y matrix nonzero. This does not affect the above argument so 
that the KK decomposition confirms the observation that the amphtude with a bulk Higgs is also 



finite [H]. Of course, for a bulk Higgs the power counting in Section 6.2 gives a more direct check 
of finiteness. Next, note that without arguing the nature of the zeros in the gauge basis Yukawa 
matrix or the physical nature of the e mixing with KK modes, it may appear that the 1/Mkk 
dependence of yijyj2 requires a "miraculous" fine tuning between the matrix elements of ( |6.11[ ). 
Our discussion highlights the physical nature of this cancellation as the mixing with heavy states 
that is unaffected by the 0{1) mixing of light chiral states. 

Finally, let us point out that the above arguments are valid for the neutral Higgs diagram 
where y = pe, the charged lepton Yukawa matrix. The analogous charged Higgs diagram contains 
neutrino Yukawa matrices yN so that there is no additional l/k from mixing. 



6.6 Matching KK and loop cutoffs 



There is one particularly delicate point in the single-mass-insertion neutral Higgs loop in the KK 
reduction that is worth pointing out because it highlights the relation between the KK scales 
M^*^ and the 5D loop momentum. To go from the 5D to the 4D formalism we replace our 
position/momentum space propagators with a sum of Kaluza-Klein propagators, 



N 



ra=0 



(6.13) 



The full 5D propagator is exactly reproduced by summing the infinite tower of states, N ^ oo. 
More practically, the 5D propagator with characteristic momentum scale k is well-approximated 
by at least summing up to modes with mass M„ ^ k. Modes that are much heavier than this 
decouple and do not give an appreciable contribution. Thus, when calculating low-energy, tree- 
level observables in 5D theories, it is sufficient to consider only the effect of the first few KK 
modes. On the other hand, this means that one must be careful in loop diagrams where internal 
lines probe the UV structure of the theory. In particular, significant contributions from internal 
propagators near the threshold M„ ^ k would be missed if one sums only to a finite KK number 
while taking the loop integral to infinity. This is again a concrete manifestation of the remarks 



below (6.1) that the length scales probed by a process depend on the characteristic momentum 



scale of the process. 

Indeed, a Kaluza-Klein decomposition for a single neutral Higgs yields 



\M\ 



ip+p')" 



Co + Cl 



M 



+ 



v y 

m) 



(6.14) 
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for some characteristic KK scale M ^ Mkk and dimensionless coefficients Cj that include a loop 
integral and KK sums. In order to match the 5D calculation detailed above, we shall work in the 
mass insertion approximation so that there are now two KK sums in each coefficient. The leading 
Co term is especially sensitive to the internal loop momentum cutoff A relative to the internal KK 
masses, 

n=l m=l \ ' J \ 'l n=\ m=l 

where we have written mass scales in terms of dimensionless numbers with respect to the mass 
of the first KK mode: M„ ~ hMkk and A ~ AMkk- It is instructive to consider the limiting 
behavior of each term c{n, m) for different ratios of the KK scale (assume n = m) to the cutoff 
scale A: 

—j for n ^ A (6.16) 

/"fl\0 

co(n, n) — T' ( — ) for n ^ A (6-17) 



.A 

Coin.n) — > ( — J for X. (6.18 



n 



We see that the dominant contribution comes from modes whose KK scale is near the loop mo- 
mentum cutoff while the other modes are suppressed by powers of the ratio of scales. In particular, 
if one calculates the loop for any internal mode of finite KK number while taking the loop cutoff 
to infinity, then the Cq contribution vanishes because the n ~ A contributions are dropped. From 
this one would incorrectly conclude that the leading order term is Ci and that the amplitude is 
orders of magnitude smaller than our 5D calculation. Thus one cannot consistently take the 4D 
momentum to infinity without simultaneously taking the 5D momentum (i. e. KK number) to 
infinity. Or, in other words, one must always be careful to include the nonzero contribution from 



modes with n ^ X. One can see from power counting on the right-hand side of (6.15) that so long 
as the highest KK number and the dimensionless loop cutoff A are matched, Co gives a nonzero 
contribution even in the A — )■ oo limit. 

This might seem to suggest U V sensitivity or a nondecoupling effect. However, we have already 
shown that /i — )■ 67 is UV-finite in 5D. Indeed, our previous arguments about UV finiteness tell 
us that the overall contribution to the amplitude from large loop momenta (and hence high KK 



numbers) must become negligible; we see this explicitly in the UV limit of (6.15). The key 
statement is that the KK scale and the UV cutoff of the loop integral must be matched, N > X. 
This can be understood as maintaining momentum-space rotational invariance in the microscopic 
limit of the effective theory (much smaller than the curvature scale). Further, the prescription that 
one must match our KK and loop cutoffs A^ > A is simply the statement that we must include all 
the available modes of our effective theory. It does not mean that one must sum a large number 
of modes in an effective KK theory. In particular, one is free to perform the loop integrals with a 
low cutoff A ~ Mkk so that only a single KK mode runs in the loop. This result gives a nonzero 
value for co which matches the order of magnitude of the full 5D calculation and hence confirms 
the decoupling of heavy modes. 
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Figure 10: Yin-Yang and double rainbow topologies of two-loop diagrams. The dotted line repre- 
sents either a gauge or Higgs boson. We have omitted the photon emission and an odd number 
of mass insertions. 



6.7 Two-loop structure 

As with any 5D effective theory, the RS framework is not UV complete. This nonrenormahzabihty 
means that it is possible for processes to be cutoff-sensitive. Since an effective /x — >• 67 operator 
(in the sense of Appendix |A]) cannot be written at tree level, there can be no tree-level counter 
term and so we expect the process to be finite at one-loop order, as we have indeed confirmed 
above. In principle, however, higher loops need not be finite. 

The one-loop analysis presented thus far assumes that we may work in a regime where the 
relevant couplings are perturbative. In other words, we have assumed that higher-loop diagrams 
are negligible due to an additional (yf^/lGvr^ suppression, where g is a generic internal coupling. This 
naturally depends on the divergence structure of the higher-loop diagrams. If such diagrams are 
power-law divergent then it is possible to lose this window of perturbativity even for relatively low 
U V cutofF A ~ Mkk- We have shown that even though naive dimensional analysis suggests that 
the /i — 67 amplitude should be linearly divergent in 5D, the one-loop amplitudes are manifestly 
finite. 

Here we argue that the two-loop diagrams should be no more than logarithmically divergent 
for bulk bosons so that there is an appreciable region of parameter space where the process is 
indeed perturbative and the one-loop analysis can be trusted. This case is also addressed in 



The relevant topologies are shown in Fig. 10 In this case, the power counting arguments that we 



have developed in this section carry over directly to the two-loop diagrams: 

Loop integrals (d^k) +8 

Gauge invariance {p + p') —1 

Bulk boson propagators —2 

Bulk vertices (dz) —5 



Total degree of divergence O 

We find that the superficial degree of divergence is zero so that the process is, at worst, logarith- 
mically divergent. 

The power counting for the brane-localized fields is more subtle, as we saw above. Naive power 
counting suggests that the two-loop, brane-localized diagrams are no more than quadratically 
divergent. However, just as additional cancellations manifested themselves in the one-loop, brane- 
localized case, it may not be unreasonable to expect that those cancellations might carry over to 
the two-loop diagrams. Checking the existence of such cancellations requires much more work we 
leave this to a fuU two-loop calculation. 
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7 Outlook and Conclusion 



We have presented a detailed calculation of the jj, ^ e^y amplitude in a warped RS model using 
the mixed p osit ion/momentum representation of 5D propagators and the mass insertion approx- 
imation, where we have assumed that the locahzed Higgs VEV is much smaller than the KK 
masses in the theory. Our calculation reveals a new sensitivity to the specific flavor structure of 
the anarchic Yukawa matrices since this affects the relative signs of coefficients that may interfere 
constructively or destructively. We thus find that while generic flavor bounds can be placed on 
the lepton sector of RS models, one can systematically adjust the structure of the Ye and 
matrices to evade the bounds while simultaneously maintaining anarchy. In other words, there 
are wide regions of parameter space which satisfy experimental constraints without fine tuning. 
Conversely, it is easy to generate anarchic flavor structures which — for a given KK scale — cannot 
satisfy the /i — )• 67 constraints for any value of the anarchic scale Y^,. Over a range of randomly 
generated anarchic matrices, the parameter controlling this Kf-independent structure has a mean 
value of zero and a la value which can push the KK scale to 30 TeV. 

It is interesting to consider the case where Mkk = 3 TeV where KK excitations are accessible 
to the LHC When the b coefficient takes its statistical mean value, 6 = 0, the minimal model 
suffers a O {10) tension between the tree- level lower bound on F* and the loop- level upper bound, 

r* > 4 Y,< 0.5. (7.1) 

This tension is slightly alleviated in the custodial model, 

K > 1.25 K < 0.3. (7.2) 

Thus for Mkk = 3 TeV one must one must accept some mild tuning in the relative sizes of the 
5D Yukawa matrix. 

On the other hand, we know that anarchic models generically lead to small mixing angles (see 
however [I7j). These fit the observed quark mixing angles well but are in stark contrast with the 
lepton sector where neutrino mixing angles are large, 0{1), and point to additional flavor structure 
in the lepton sector. For example in [15] a bulk A4 non-Abelian discrete symmetry is imposed on 
the lepton sector. This leads to a successful explanation of both the lepton mass hierarchy and 
the neutrino mixing angles (see also [32J) while all tree- level lepton number-violating couplings 
are absent, so the only bound comes from the /i — )■ 67 amplitude. 

We have also provided different arguments for the one-loop finiteness of this amplitude which 
we verified explicitly through calculations. We have illuminated how to correctly perform the 
power counting to determine the degree of divergence from both the 5D and 4D formalisms. The 
transition between these two pictures is instructive and we have demonstrated the importance of 
matching the number of KK modes in a 4D EFT to any 4D momentum cutoff in loop diagrams. 
The power-counting analysis can be particularly subtle for the case of brane-localized flelds and we 
have shown how one-loop finiteness can be made manifest. Finally, we have addressed the existence 
of a perturbative regime in which these one-loop results give the leading result by arguing that 
the bulk field two-loop diagrams should be at most logarithmically divergent and that it is at least 
feasible that the brane-localized two-loop diagrams may follow this power counting. 

In addition to fi e'j is an analogous process mediated by the same penguin diagrams, b — )• 57. 
There are additional gluon diagrams with the same topology as the Z diagrams described here. 
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Because of operator mixing, connecting the 6 — )• 57 amplitude to QCD observables requires the 
Wilson coefficients for both the photon penguin C?^ and the gluon penguin Csg. A discussion 
can be found in though there it was expected that these penguins would be logarithmically 
divergent. We leave the exphcit evaluation of the 6 — )■ 57 amphtude in warped space to future 
work [33] . 
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A Matching 5D amplitudes to 4D EFTs 

The Standard procedure for comparing the loop-level effects of new physics on low-energy ob- 
servables is to work with a low-energy effective field theory in which the UV physics contributes 
to the Wilson coefficient of an appropriate local effective operator by matching the amplitudes 
of fuU and effective theories. In this appendix we briefly remark on the matching of 5D mixed 
position/momentum space amplitudes to 4D effective field theories, where some subtleties arise 
from notions of locality in the extra dimension. 

The only requirement on the 5D amplitudes that must match to the 4D effective operator is 
that they are local in the four Minkowski directions. There is no requirement that the operators 
should be local in the fifth dimension since this dimension is integrated over to obtain the 4D 
operator. Thus the 5D amplitude should be calculated with independent external field positions 
in the extra dimension. Heuristically, one can write this amplitude as a nonlocal 5D operator 

05(0;, zh, zl, ze, za) = H^{x, zh) ■ L^{x, zl) o"^^^ E^{x, ze)Fmn{x, za). (A.l) 

Note that this object has mass dimension 8. In the 5D amplitude the fields are replaced by external 
state wavefunctions and this is multiplied by a "nonlocal coefficient" c^{zh,zl, ze, za) which in- 
cludes integrals over internal vertices and loop momenta as well as the mixed position/momentum 
space propagators to the external legs. To match with the low-energy 4D operator we impose that 
the external states are zero modes and decompose them into 4D zero-mode fields multiplied by a 
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5D profile f (z) of mass dimension 1/2, 

$5(X,2) ^ <l>(°)(x)/(°)(2). 



(A.2) 



Further, we must integrate over each extemal field's z-position. Thus the 4D Wilson coefficient 
and operator are given by 



Wdzi 



C,{zH.ZE.ZL,ZA)fP{zH)f'E\zE)fr{zL)fT{zA) H ■ La^^ EF,,, (A.3) 



p(0). 



(0), 



where the fields on the right-hand side are all zero modes evaluated at the local 4D point x. Note 
that these indeed have the correct 4D mass dimensions, [C4] = 6 and [c] = —2. 

Finally, let us remark that we have treated the 5D profiles completely generally. In particular, 
there are no ambiguities associated with whether the Higgs field propagates in the bulk or is 
confined to the brane. One can take the Higgs profile to be brane-locahzed, 



fn^zn) ~ ^^'5{z - R') 



(A.4) 



where the prefactor is required by the dimension of the profiles. With such a profile (or any 
limiting form thereof) the passage from 5D to 4D according to the procedure above gives the 
correct matching for brane-localized fields. 



B Estimating the size of each diagram 

As depicted in Figs. |2]-[4| there are a large number of diagrams contributing to the a and b coef- 
ficients even when only considering the leading terms in a mass-insertion expansion. Fortunately, 
many of these diagrams are naturally suppressed and the dominant contribution to each coeffi- 
cient is given by the two diagrams shown in Fig. [5j This can be verified explicitly by using the 
analytic expressions for the leading and next-to-leading diagrams are given in Appendix [C} In this 
appendix we provide some heuristic guidelines for estimating the relative sizes of these diagrams. 



B.l Relative sizes of couplings 



First note that after factoring out terms in the effective operator in (4.2), Yukawa couplings give 
order one contributions while gauge couplings give an enhancement of In R'/ R, where gsu is 
the appropriate Standard Model coupling. This gives a factor of ~ 5 (7) enhancement in diagrams 
with a W over those with a Z (H). 



B. 2 Suppression mechanisms in diagrams 

Next one can count estimate suppressions to each diagram coming from the following factors 

A. Mass insertion, ~ 10^^/insertion. Each fermion mass insertion on an internal line intro- 
duces a factor of O {v R'). This comes from the combination of dimensionful factors in the 
Yukawa interaction and the additional fermion propagator. 
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Bi. Equation of motion, ~ 10"'^. Higgs diagrams without an explicit chirality-flipping internal 
mass insertion must swap chirality using the muon equation of motion u{p)p = m^u{p). This 
gives a factor of 0{m^R') and is equivalent to external mass insertion that picks up the zero- 
mode mass. 

B2. External mass insertion, ~ 10~^. Alternately, when a loop vertex is in the bulk, an exter- 



nal mass insertion can pick up the diagonal piece of the propagator — see (G.l ) — representing 
the propagation of a zero mode into a 'wrong-chirality' KK mode. Unlike the off-diagonal 
piece which imposes the equation of motion, this is only suppressed by the O {v R') mentioned 
abov^ One can equivalently think of this as an insertion of the KK mass which mixes the 
physical zero and KK modes. 

C. Higgs/Goldstone cancellation, ~ 10^'^. The and one-mass-insertion loops cancel 
up to O {{jn^j — 'rn?z) / ml^Y^ because the two Goldstone couplings appear with factors of i 
relative to the neutral Higgs coupling^ 

D. Proportional to charged scalar mass, ~ 10^^. The leading loop-momentum term in 
the one-mass-insertion brane-localized loop cancels due to the form of the photon cou- 
pling relative to the propagators. The gauge-invariant contribution from such a diagram is 



proportional to {My/B!) . This is shown explicitly in (B. 2) below 



To demonstrate the charged scalar mass proportionality, we note that the amplitude for the 
one mass insertion charged Higgs diagram in Fig. |3] is 

M^--E'(^\'^f Y'f f v A^A^v {2k-p-pr 

(B.l) 

Remembering that the 5D fermion propagators go like A ~ ^/k, this amplitude naively appears to 
be logarithmically divergent. However, the Ward identity forces the form of the photon coupling 
to the charged Higgs to be such that the leading order term in fc^ cancels. This can be made 
manifest by expanding the charged Higgs terms in p and p', 



{2k-p- p'Y {p + p') 



P _ p/)2 _ M^] [{k - p)^ - M^] (P - M, 



2 

W) 



1 



where we have dropped terms of order (9(m^/M^). Thus see that the coefficient of the gauge- 
invariant contribution is finite by power counting. After Wick rotation, this amplitude takes the 
form 

M^{IM\H^)\^^^^,^ = —{B!)^f,^Y^f_,^^ — {E!Mw)'hMm^u,,{p + p')u,, (B.3) 



where /imih± is a dimensionless integral given in (jCj). We see that the amplitude indeed carries a 
factor of [MwR!)'^. 



■^We thank Martin Beneke, Paramita Dey, and Jiirgen Rohrwild for pointing this out. 
^We thank Yuko Hori and Takemichi Okui for pointing this out. 
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Figure 11: Alignment of the extemal mass insertion diagrams with Standard Model gauge bosons. 
X and ip are left- and right-chiral Weyl spinors respectively. The gauge boson vertices don't 
change fermion chirality so that the internal fermion must be a chirality-fiipping KK mode. We 
have neglected the contribution where the external mass insertion connects two zero mode fermions 
since this is suppressed by m^R' . 



B. 3 Dimensionless integrals 

Estimating the size of dimensionless integrals over the loop momentum and bulk field propagators 
(such as /iMiH±) is more subtle and is best checked through explicit calculation. However, one 
may develop an intuition for the relative size of these integrals. 

Note that the fifth component of a bulk gauge field naturally has boundary conditions opposite 
that of the four-vector [20] so that the fifth components of Standard Model gauge fields have 
Dirichlet boundary conditions. This means that diagrams with a W^H^A vertex vanish since 
the brane-localized Higgs and bulk do not have overlapping profiles. Further, loops with fifth 
components of Standard Model gauge fields and internal mass insertions tend to be suppressed 
since the mass insertions attach the loop to the IR brane. In the UV limit the loop shrinks towards 
the brane and has reduced overlap with the fifth component gauge field. 

Otherwise the loop integrals are typically C(O.l). The particular value depends on the propa- 
gators and couplings in the integrand. 



B. 4 Robustness against alignment 



As discussed in Section 5.2, the flavor structure of the diagrams contributing to the h coefficient 
is aligned with the fermion zero-mode mass matrix |H [TU |T2] . Contributions to this coefficient 
vanish in the zero mode mass basis in the absence of additional fiavor structure from the bulk 
mass (c) dependence of the internal fermion propagators. The diagrams which generally give the 
largest contribution after passing to the zero mode mass basis are those with with the strongest 
dependence on the fermion bulk masses. Since zero mode fermion profiles are exponentially 
dependent on the bulk mass parameter, a simple way to identify potential leading diagrams is to 
identify those which may have zero mode fermions propagating in the loop. 

This allows us to neglect diagrams with an external mass insertion and a 4D vector boson in 



the loop. As shown in Fig. 11, such diagrams do not permit intermediate zero modes to leading 



order. Note, however, that diagrams with an external mass insertion and the fifth component of 



gauge boson are allowed to have zero mode fermions in the loop. Indeed, in Section 5.2 we show 
that such a diagram (with a and W^) is the dominant contribution in the minimal model since 
other diagrams {Z with an internal mass insertion, H^/W^) end up being numerically small. 
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c Analytic expressions 



We present analytic expressions for the leading and next-to-leading diagrams contributing to 
/i — )■ 67. We label the diagrams in Figs. [2]-[4] according to the number of Higgs-induced mass 
insertions and the internal boson(s). For example, the two-mass-insertion W diagram in Fig.|5|i 
is referred to as 2MIVr. Estimates for the size of each contribution are given in Appendix B] 
We shall only write the coefficient of the Upr{p + p')^Up term since this completely determines the 
gauge-invariant contribution. 



Cl Dominant diagrams 

As discussed in Section |5| the leading diagrams contributing to the a and b coefficients are 
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We have explicitly labeled the 4D (dimensionless) anarchic Yukawa matrices whose elements as- 
sumed to take values of order {YE)ij ~ (^Af)ij ~ Y*i but have independent flavor structure. Note 
that we have suppressed the flavor indices of the Yukawas and the dimensionless integrals. Dia- 
grams with a neutral boson and a Yukawa structure YeY^Ye also contribute to the a coefficient, 
but these contributions are suppressed relative to the dominant charged boson diagrams above. 
These diagrams may become appreciable if one permits a hierarchy in the relative Ye and Y^ an- 
archic scales, in which case one should also consider the Z boson diagrams whose analytic forms 
are given below. 

The dimensionless integrals are 
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where x = kEZ, y = kER', and = m^R! . The significance of these dimensionless variables is 
discussed below (6.1). The dimensionless Euclidean-space propagator functions F are defined in 



(F. 39 - F. 40), where the upper indices of the F functions define the propagation positions. For 
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example, F^^"^ represents a propagator from z = R' to z = z^. Similarly, G y and G y are defined 



in (E.4) and (E.5) 



C. 2 Subdominant a coefRcient diagrams 

The diagrams containing a brane-localized Higgs loop are 



M{nMlH^) 



{p+p')" l67r 
i 



MI_ff±! 



ev 



(C.5) 
(C.6) 



Here n = 0, 1 counts the number of internal mass insertions in the diagram. The gauge boson 
loops are 



M{nMlZ/Z^^ 



(P+P'Y ~ I67r2 fcLYsYEYEf-CE ^ 



(^^n'lnf (^i?')'x W/z., 

(C.7) 



Mi2Mlww)\^p_ 



Where n = 2, (1 + 2), 3 with (1 + 2) referring to a single internal mass insertion and two external 
mass insertions. 2M\ww represents 2MW^iy^, 2MWiy^ and 2MIiy^iy. The dimensionless 
integrals are 



/ij^i^o = / dydxy^ f'^V f-2FfrFf?F^^^ 



+,y +,y -,y „,2 



2/2 + {MhR')^ 



I pLyx pLxy pB.yy ^ 1 /^7, r) P^y^\ pLxy pB.yy 

+ ^+,, ^+,, J'-,y + (rnnE')^)^ 2 'v^'^'^^^+'f J ^+'?^ ^"'S' 



2/2 + (MhR') 



- f v 3, n F^y'^^ f) . F^'^y F^y^ ^ \-2F^yy n . F^y^" n F^^'y 3^ 



_ pLyy jj pRyx jj pRxy y_ 

+ 'y (2/2 + {MhR'W ' 2 

pLyy pRyx pRxy ^ |^ 1 TP^yy\ fpRy^ fpR^y ^ 

1 ^ 



\ _ { ^.f) p^yy \ ptiyx pnx 

y --,y - -.y ^2 + [MuRIf 2 V'^'^^^^^y J ^"-^^ ^--J/ y2 + (M^i?')' 



.1! L^F^yy (vd, D.F^y'^D F^'^y — 

?A2 + o^+.y ] „2 ^ rA#..R/^2^ 



2/2 + {M H R' f 2 



2/2 + (Mh/2')' 
(C.9) 
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-^lMI_ff± 



'2MIZ'' 



'(H-2)MIZ 



2 y' 



Ryy 



-I 



huiiz ^ dy dxidx2dx3 



XiJ 
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2+cl / \ 4 / \ 4 



y_ 

X2 



X 



{y d,, Gf D^Fily; [y' [f^^ f^' F^^y d.f':^' + F'^y fZ' D.F!^f 



I pLly fpRyy f) fpR2y pL23 , pLly fpRyy fpLy2 r^L23\ 

'^^+,y ^-,y ^-^-,y ^-,y ^ +,y ^ -,y ^+,y ^-^-,y J 



-,y 



-,y 



D F^^'^ D^F^'^y F^yy d P^y^ 

^-^-y ^+^+,y ^-,y ^-^-,y 



-( 

— j dy dxidx2dx3 

^ y d,^ Gf D^Ftly; (y' - ^^F^f/ F^-) 
dy dxidx2dxy, 



Fi:\y D^F^f b.F'^^y D+P'-y^ 



y_ 

X2 



,y 

y_ 

X3 
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X 



y_ 

Xi 



y_ 

X2 



D^p;i:' D^Pily; Gf - (4 + ^ 9. J [D.P^^y P^fn.P^^f + Ff/ D_F 
+ D.P'^ly F»_Ff,f F^f ^ f)_pR^y pLys f)_pL^^ 



R3y 
y 



-^(i+2)Miz^ ~ ~ / dy dxidx2dx 



y_ 
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y_ 

X2 



D^PZf D^P^ly; Gf 
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PZ' D_PZ y d,, (i)^F-3) + , 4, (y' P^]^ P^) P^ 
The integral for 3MIZ and SMIZ^ can be written as 

2+c_c / _ N 2-CE 



'3MIZ/Z5 



dy dxidx2dx3 



y_ 

Xi 



y_ 

X3 



Gf^M,ydk,Ni. 
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For 3MIZ, the (M, N) pairs are 

M2 
M3 
M4 

Mg 

Mg 
M7 



A^2 

N, 

N7 

N. 



For 3MIZ^ the (M, A^) pairs are 



iVa) = 

A^s) = 
N7) = 



pL12 4 pL2y pRyy fpLyy pRy3\ 

^+,y ' y ^+,y ^-,y ^+,y ^-,y J ' 

-y' D+Ft^ F^l' Fil' F^f , D-F^l]^) , 

_„,2 pR2y ^Lyy ^Ry3 _ 2 ^Lly 

y ^~,y ^+,y ^-,y ; y ^+,y ^~,y ^ 

^+,y ^+^+,y ' y ^-^-,y ^+,v ^ -,y ) ' 

—1/2 W^'^y P^yy P^y"^ p-^^y ^RyS^ 

y ^ +,y ^-,y ^ +,y 1 y ^ +,y -,y 

f) , pL2y pRy3 2 pLly pRyy f, pLy2\ 

pR23 „,4 pLly fpRyy fpLyy pRy2\ 

^-,y ' y ^+,y ^-,y ^+,y ^-,y J ' 
y ^+,y ^-,y ^+,y ^+^+,y ' ^-^-,y J ■ 

-y' b^Fily P^^y Piiy P^f , D^P^f) , 

FZ' ' -y' D+F'c'y P^'y D+F^^) , 

b^Fi^y F^y b_F!:i\ b^Fil^ d^fZ') 

p^^b^p^f , -y' b^Fii' fZ' Fii') , 

D-F^^ fZ' D^fZ' , b^F^'D^FZ') 

F^, -y' D^F^^ fZ' F^y D+F^;; 



The integrals for the loops are 
-^2Miw^w^ ~ — I dy dxidx2dx2, 



-7/2 P^'^y P^yy P^yy b , P^y^ b , P^'^A 

y ^+,y ^-,y ^+,y ^+^+,y ' ^+^+,y ) ■ 



y_ 

X2 



y_ 

X3 



X 



1 



y^ D^Ft:,' F-yy FZ' D^F-y; 4 G- G- + yd,l G^' 



-^2MIW5W 



-/ 



dxidx'2,dx^ 



Xi 
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X 



(C.34) 



y2 pLly pRyy pLy2 p^^pL2y^ ^13 Q^ ^23 _ y ^13 ^23^ j (^.35) 
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-^2Miww5 = - / dydxidx2dx3 [ — j I— j ( — ) x 



y \ I y\ y 



2 



D^Fil^ F^^y D.F!:f D^F^l^ [y d,, Gf d. Gf - y d,,d. Gf Gf) 



C. 3 Subdominant h coefRcient diagrams 

M (OM W) I (^^^,), = ^ (R'f f,^YEf-c,^ {gl^)' In^ X /omw, 
where n = 0, 1 counts the number of internal mass insertions. 

// \ 2+C£, / \ 2-ce / \ 4 
dydxidx2dx3 i —j i —j i —j {y dk^G^^) y'^ x 

( - D^F^f D.F^^y Fif + F^J- pR^y pLyi 



_ pR2y f. pLy3 f\ pL31 , pR2y pLy3 fpLZl 
^-,y ^-^-,y ^+^+,y ~^ ^-,v ^+,y ^+,y 

2+CL / \ 2-CB / \ 4 
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\ [Ff y [D^p^y D^F^f) Gf - D^F^^ y d,^ (f^ J D^F^ G 

- 4 D^F^]^ F^'y b^F^f Gf + D+F^ J F^f (y d,, D^F^f) Gf 

- b^Fi^y b^PZ' y 4. (f^' Gf) - 4 b^Fily b^PZ' P^f g f 

-^oMiz = dy dxidx2dx3 ( — ) ( — ) ( — ) x 

y d,^ G f b^P^y^ (^_Ff - F-3 + ^.12 ^_^.23) 

ToMiz^ ^ - J dy dxidx2dx3 x 

{ i " + y d,, Gf) + y d,, Phf G'f\ } , 

loMiw^- J dydxidx2dxs ^-^^ ^-^^ ^-^^ x 
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C.4 Custodial Models 

For custodially protected models, one must include loops with the custodial partners of fermions 
and gauge bosons. See, e. g., [S| for details of the additional field content of such models. The new 
particles have mixed boundary conditions, ( — h) or (H — ). For the chirality flipping process fi — )■ 
6 7, Yukawa insertions on the IR brane only allow fermions carrying either (++) or ( — h) boundary 
conditions running in the loop. This limits the number of the new diagrams to be considered. 
The new fermion propagators can be obtained by making the replacement F ^ E. Writing the 
boundary condition in terms of the Weyl components of the Dirac spinor, corresponds to 
the boundary condition (■?/'(+_), , while E^ corresponds to (■?/'(_+), • For x > x', the 

£'-functions can be written as foUows: 



{xx'f/^Sc{x.,y^)T^{x_,wy+) {xx'fl^T^{x+,y.)S^{x^,wy+) 

r Tc{y-,wy+) y'> Tc[y-,wy+) 

{xx'f/^T^{x_,y+)Sc{x'_,wy_) {xx'f/^ Sc{x+,y+)T^{x'^,wy_) 

r Tc{y+,wy-) y^ Tc{y+,wy^) 

The X < x' expressions are obtained by replacing x -H- x'. Gauge bosons with ( — h) boundary 
conditions can also appear in custodial loops. The corresponding propagator for x > a;' is G — )• iif 
with 

n (Rf xx' T,o{x,y)Suix',wy) 

hk[x,x) = — — — r . (L.48j 

R y Tio{wy,y) 

The T and S are defined in Appendix. ([E|), and the x < x' case can be obtained by x -H- x'. 



D Position, momentum, and position/momentum space 

In order to elucidate the power counting in Section [6] and to provide some motivation for the 
structure of the propagators in Appendix |F.1[ we review the passage between Feynman rules in 
position, momentum, and mixed position/momentum space. For simplicity we shall work with 
massless scalar fields on a flat (Minkowski) ci-dimensional background, but the generalization of 
the salient features to higher spins is straightforward. In position space, the two-point Green's 
function for a particle propagating from x' to x is 

D{x,x') = j d'^k^e~'^-^''-'''\ (D.l) 

a momentum-space integral over a power-law in k times a product of exponentials mk-x and k-x' . 
Each vertex carries a df'x integral representing each spacetime point at which the interaction may 
occur. When some dimensions are compact, the associated integrals are reverted to discrete sums 
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and the particular linear combination of exponentials is shifted to maintain boundary conditions. 
Further, when dimensions are warped the exponentials become Bessel functions. In this Appendix 
we will neglect these differences and focus on general features since the UV behavior of each of 
the aforementioned scenarios (i. e. for momenta much larger than any mass, compactification, or 
warping scales) reduces to the flat noncompact case presented here. 

In 4D it is conventional to work in full momentum space where the Feynman rules are derived 
by performing the df'x integrals at each vertex over the exponential functions from each propagator 
attached to the vertex and amputating the external propagators. This generates a momentum- 
conserving (5-function at each vertex which can be used to simplify the d'^k integrals in each 
propagator. For each diagram one such 5-function imposes overall conservation of the external 
momenta and hence has no dependence on any internal momenta. For a loop diagram this means 
that there is a leftover d'^k which corresponds to the integration over the loop momentum. Thus 
the momentum space formalism involves separating the exponentials 'm k ■ x from the rest of the 
Green's function and performing the d'^x integral to obtain 5-functions. 

To go to the mixed position/momentum space formalism we pick one direction, z, and leave 
the dependence on that position in the propagator while integrating over the z-component of the 



momentum, k^ in (D.l). We shall write the Minkowski scalar product of the {d — 1) momentum- 



space directions as k so that the full (i-dimensional scalar product is k — k^. The Feynman rule 
for each vertex now includes an explicit dz integral which must be performed after including each 
of the position/momentum space propagators, which take the form 

A(fc,z,/) = j dk,^-^^e'''^^'-''\ (D.2) 

The {d — 1) other exponentials and momentum integrals are accounted in the usual momentum- 
space formalism. This object goes like A ~ which indeed has the correct dimensionality for 
the sum over a KK tower of scalar propagators. Similarly, the massless bulk fermion propagator 
is 

A(A:, z, z') = f dh^f^e^'^^-^l (D.3) 



where we may now identify the scalar functions F ~ (i/c^e*^^*-^^^''*/(^^ — k^) in (F. 6) and (F. 23). 

It is thus apparent that the mixed formalism contains all of the same integrals and factors as 
the momentum-space formalism, but that these are packaged differently between vertex and prop- 
agator Feynman rules. By identifying features between the two pictures one may glean physical 
intuition in one picture that is not manifest in the other. For example, the observation in the 
mixed formalism that each bulk vertex on a loop brings down a power of l/fc is straightforwardly 
understood to be a manifestation of momentum conservation in the momentum space picture. 

On the other hand, the mixed formalism is much more intuitive for brane-localized effects. 
Interactions with fields on the brane at z = L carry 6{z — L) factors in the vertex Feynman rules. 
Such interactions violate momentum conservation in the 2;-direction. In the KK formalism this 
manifests itself as the question of when it is appropriate to sum over an independent tower of 
KK modes. This is easily quantified in the mixed formalism since the dz integrals are not yet 
performed in the Feynman rules and we may directly insert 6{z — L) terms in the expression for 
the amplitude. 
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Figure 12: A simple loop diagram to demonstrate the power counting principles presented. The 
hnes labeled pi represent the net external momentum fiowing into each vertex so that pf corre- 
sponds to the KK mass of the i^^ external particie. 



As a concrete example, consider the loop diagram with three vertices shown in Fig. [12} It is 
instructive to explicitly work out loop z-momentum structure of this diagram in the case where 
all vertices are in the bulk and observe how this changes as vertices are localized on the brane. 
To simplify the structure, let us define the product of momentum-space propagators 

3 

f{h,h,h)^l[ ^,' (D.4) 

i=l * \ i I 

Using J dz exp{izk) = 6{k), the bulk amplitude is proportional to 

j dzi dz2 dzs dkl dkl f{k^, fca, k^) e'^'^^'+P'-'"''^^ ^iz2ik2+P3-k,r ^iz^{k,+pz-k,r (135) 

dz2 dzs dkl dkl fik2 - Pl, ^2, ^3) e^"2(^2+P3-*^^)' ^iz^ik^+P^-k^+p^Y (d.6) 

dzs dkl f ih -P2- Pl, ks - P2, ks) e*^^(Pl+^'2+f3)^ (D.7) 



We have implicitly performed the associated d^'^^^^z integrals at each step. The final dzs integral 
gives the required 5-function of external momenta while leaving an unconstrained dkl loop integral. 
Each dk^ /{k"^ — k1) ~ l//c represents the entire KK tower associated with an internal line. The 
removal of two dk'^ integrals by (5-functions is a manifestation of the l/k suppression coming from 
each dz integral with the caveat that the "last" dz integral only brings down powers of external 
momenta and hence does not change the power of loop momenta. This explains the "overall 



2;-momentum" contribution to the superficial degree of divergence in Section |6.2 

Next consider the case when the zs vertex is brane localized so that its Feynman rule is 
proportional to 5(^3 — L). This only affects the last line of the simplification by removing the dzs 
integral. Physically this means that 2r-momentum (KK number) needn not be conserved for this 
process. Since the ^3 exponential is independent of any loop momenta, this does not affect the 
superficial degree of divergence. 

On the other hand, if Z2 is also brane localized, then the 5{z2 — L) from the vertex prevents the 
dz2 integral in the second line from giving the 5{k2 + P2 — ks) that cancels the dk^ integral. Thus 
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the process has an additional integral which now increases the degree of divergence. In the 
4D formalism this is manifested as an additional independent sum over KK states. It is now also 
clear that setting zi to be brane localized prevents the dk^ from being cancelled and hence adds 
another unit to the degree of divergence. This counting is trivially generalized to an arbitrary 
number of vertices and different types of internal propagators. For a loop with V vertices, Vb of 
which are in the bulk, the key points are: 

1. U V = Vb, then the dz integrals reduce the superficial degree of divergence by {Vb — 1). 

2. If, on the other hand, \^ > Vb so that there is at least one brane-localized vertex, then the 
dz integrals reduce the superficial degree of divergence by Vb- 

Intuitively the z-momentum nonconservation coming from brane-localized interactions can be 
understood as the particle picking up an arbitrary amount of momentum as it bounces ofF the 
brane (a similar picture can be drawn for the orbifold [34j). Alternately, it reflects the uniform 
spread in momentum associated with complete localization in z-position. While this may seem 
to imply sensitivity to arbitrarily high scale physics on the brane, a negative degree of divergence 
will prevent the loop from being sensitive to UV physics. In other words, we are free to treat 
brane-localized fields as having (5-function profiles independent of the physics that generates the 
brane. 

Finally, note that we have assumed that each fermion mass insertion is brane localized. In 
5D this means that higher-order diagrams in the fermion mass-insertion approximation are not 
suppressed by momentum since each additional brane-to-brane propagator goes like ~ ^//c after 
accounting for the dk^ integrals. Instead, these mass insertions are suppressed only by the relative 
sizes of the Higgs vev and compactification scale, {vR'^ ~ .01. It is perhaps interesting to note 
that our analysis further suggests that in 6D with a Higgs localized on a 4D subspace, there are two 
additional momentum integrals coming from a mass insertion so that each vev-to-vev propagator 
goes like a positive power of the momentum ~ ^ causing the mass-insertion approximation to 
break down. 

E Bulk Feynman Rules 

Here we summarize the 5D position/momentum space Feynman rules used to derive the amplitudes 
in this paper. AU couplings are written in terms of 5D quantities. The brane-localized Higgs field 
is drawn as a dashed line and the fifth component of a bulk gauge boson is drawn as a dotted line. 
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'i?' 



Ak(z,z') 



/VVN/V^ = -ir]'"'Gk{z, z') 




-^595 V V 



iGk(z,z') 



fc (Z 



u{p 



R' \RJ \R' 

/\ f c 



R' \RJ \R' 



The 5D Lagrangian parameters are related to the usual Standard Model parameters by 

gl = gl^RlnR'/R 

65 = esM 



(E.l) 
(E.2) 

Y, = RY, (E.3) 
where Y represents an anarchic 4D Yukawa matrix that is related to the Standard Model Yukawa 



by (2.7). The /c fermion flavor functions are defined in (2.5). The vector propagator functions 



Gk{z,z') and Gk{z,z') are explicitly derived in [35], which also contains generic formulae for 
analogous functions for fields of general spin and additional gauge boson vertices. Using the 



dimensionless x and y variables defined in (6.1) and assuming z > z', the Euclidean space vector 
Green's functions are 



Gk{z,z') 
Gk(z,z') 



{R! f _ {R'Yxx'T,o{x,y)T,o{x',wy) 

(_Tjy X j - 



R 
R 



R y Soo{wy,y) 
{R' f xx' 5'oo(x, y)5'oo(x^ wy) 
R y Sooiwy,y) 



(E.4) 
(E.5) 



where 



T,jix,y) = h{x)K,{y) + I,{y)Ki{x) (E.6) 

S,,{x,y) = Ii{x)K,{y) - Ijiy)K,{x) (E.7) 

and w = R/ R'. For z < z' the above formula is modified by x ^ x'. The three gauge boson 
couplings are given by 



Wf 




R 



265- [{k - k^Yr]"'' + {k- - ky?]"" + {k^ - fc-) V"] 



W: 
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z 



w: 




w; 



Here we have used the convention where all momenta are labeled by the charge of the particle and 
are flowing into the vertex. The A^W^W~ vertex is given by e^{R/ z)ri^'^[d'^ — dz)- The Euchdan 



space fermion propagator Ak{z, z') is given in (F. 41). 



F Derivation of fermion propagators 

General formulae for the scalar function associated with bulk propagators of arbitrary-spin fields 
in RS can be found in [35]. The special case of bulk fermion propagators with endpoints on 
the UV brane is presented in [36] . The Green's function equation for the general RS fermion 
propagator can be solved directly from the Strum-Liouville equation, though this can obscure 
some of the intuition of the results. Here we provide a pedagogical derivation of the 5D bulk 
fermion propagator in a flat and warped interval extra dimension. See also the discussion in 



Appendix D which relates this construction to the usual pure momentum space formalism. 



F.l Flat 5D fermion propagator 

First we derive the chiral fermion propagator in a flat interval extra dimension z G (O, L) as a 



model calculation for the warped fermion propagator which is presented in Appendix |F.2[ A 
complete set of propagators for a flat 5D interval was derived in ^TBj using finite temperature field 
theory techniques. 

We derive these results by directly solving the Green's function equations. The propagator 
from a given point x' to a another point x is given by the two-point Green's function of the 5D 
Dirac operator, 

VA{x,x') = {i^^dM -m) A{x,x') = t6'~^\x - x'), (F.l) 

where M runs over 5D indices. We shall treat the noncompact dimensions in momentum space 
and the finite dimension is in position space. In this formalism, the Green's function equation is 

{-f + id^l^ - m) A(p, z, z') = i6{z - z'), (F.2) 

where we use 7^ = diag(il2, —H2)- 

This is a first-order differential equation with nontrivial Dirac structure. To solve this equation 
we define a pseudo-conjugate Dirac operator (which is neither a complex nor Hermitian conjugate), 

V = i-i^'dM + m. (F.3) 
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Using this to "square" the Dirac operator, we can swap the Dirac equation for a simpler Klein- 
Gordon equation that is second order and diagonal on the space of Weyl spinors, 

VV^i^'-^'-' ^,_^,_,). (F.4, 

It is straightforward to solve for the Green's functions z') of the VD* operator in mixed 

position/momentum space, 



From these we can trivially construct a solution for the Green's function of (F.l), 

We solve this by separating F±{z) into pieces 

F±{p,z,z) = \ n -r ^ , 

I z, z) \i z> z 

and then solving the homogeneous Klein-Gordon equations for each and F^ . The general 
solution is 

F<'>(p,z,z') = A<'>cos(xp^) + 5<'>sin(xp^), (F.8) 

where the eight coefficients A^'^ and i?^'^ are determined by the boundary conditions at O, L and 
z' . The factor Xp is the magnitude of ps and is defined by 

Xp = a/p^ - w?. (F.9) 

We impose matching boundary conditions at z = 2;'. By integrating the Green's function 
equation (F.5) over a sliver z ^\z' — e, z' + e] we obtain the conditions 

d,F>{z')-d,F<{z')=i, (F.IO) 
F>{z')-F<{z') = Q. (F.ll) 

These are a total of four equations. The remaining four equations imposed at the branes impose 
the chirality of the fermion zero mode and are equivalent to treating the interval as an orbifold. 
We denote the propagator for the 5D fermion with a left-chiral (right-chiral) zero mode by 
(A^). We impose that the Green's function vanishes if a "wrong-chirality" state propagates to 
either brane, 

Pr A^(p,z,/)|^^^^^ = PnV F\p,z,z')\^^^^^ = O, (F. 12) 

Pl A\p,z,z')\^^^^^ = PlV F^(p,^,^')L=o,l = O' 
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Sp{L 



)SpL 



O 
O 



Cp{L 



^pZ ^plj 

c p z c p T J 



O 

-Cp(L - z') 
Sp{L - z') 
O 



qR> 



~ C p z Sp J-j 
CpZ Cplj 
~ Sp z c p J-j 

^p z s p -L 



Table 1: Flat case coefficients in (F.8) upon solving with the boundary conditions (F. 10 -F. 13) 
We have used the notation CpX = cos XpX and SpX = sin XpX. 



where Pl,r = ^(1 =F ^7^) are the usual 4D chiral projection operators. Note from (F. 6) that each 



of these equations is actually a set of two boundary conditions on each brane. For example, the 
left-handed boundary conditions may be written exphcitly as 

f!:{p,z,z') 

{d, + m)F':{p,z,z') 



=0,L 



=0,L 



O, 

o, 



(F.14) 
(F.15) 



where we have used that is arbitrary. It is well-known that only one boundary condition for 
a Dirac fermion needs to be imposed in order not to overconstrain the first-order Dirac equation 
since the bulk equations of motion convert boundary conditions for x iiito boundary conditions 
for ijj [37j. In this case, however, we work with a seconc?-order Klein-Gordon equation that does 
not mix x and ip. Thus the appearance and necessity of two boundary conditions per brane for a 
chiral fermion is not surprising; we are only converting the single boundary condition on A (p, z, z') 
into two boundary conditions for F{p, z, z'). 

Solving for the coefficients A'^'^{p,z) and B^'^(p,z) for each type of fermion (left- or right- 
chiral zero modes) one finds the results in Table [T] Using trigonometric identities one may combine 
the z < z' and z > z' results to obtaiijf] 

-icosXpiL- \z - z' \) +-f^px cos Xp (L - {z + z')) 
2Xp sm XpL 

+1 and pn = —1. The fermion Green's function can then be 



where X = {L, R} with pi 



obtained trivially from (F. 6). 



Let us remark that the leading UV behavior of a brane-to-brane propagator (where the k^'j^ 
term vanishes) goes like 

Xk 



A 



(F.17) 



F. 2 Warped 5D fermion propagator 

We now derive the chiral fermion propagator in a warped interval extra dimension following the 



same strategy as Appendix F.l The Dirac operator is obtained from the variation of the Randall- 
Sundrum free fermion action, 



Srs (fermion) 



R' 

dx I dz 

R 



Z z 



(F.18) 



^This result difFers from that of [TS] by a factor of 2 since that paper treats the compactificd space as an orbifold 
over the entire rather than just an interval [O, Tri?]. 
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where c = mR and we have integrated the left-acting derivatives by parts. The Dirac operator is 
a product of the {R/zY prefactor coming from the AdS geometry and an operator V given by 



V = i-f^'^dM -i-l^ - -. (F.19) 
z z 

We would hke to find the mixed position/momentum space two-point Green's function satisfying 



{R/zyVA{p,z,z') =i5{z-z'). 



FoUowing (F. 3) we define a pseudo-conjugate Dirac operator 

.2 



%-r + - 

z z 



and 'square' D into a diagonal second-order operator, 
(VV 



VV 



o 

VV 



VV 



± 



, c^±c 



(F.20) 



(F.21) 



(F.22) 



Next we follow (F.5) and solve for the Green's function of this squared operator in mixed posi- 
tion/momentum space where — )■ — p^, 



-{R/zfVVF{p,z,z) = - 



VV 



VV 



iS{z — z'). 



(F.23) 



The solution to the Dirac Green's function equation (F.20 ) is then given by A(p, z, z') = VF{p, z, z'). 
We shall separate F{p,z,z') into solutions for the cases z > z' and z < z' following ( F.7[ ). The 
general solution to the homogeneous equation (F.23) with z ^ z' is 



F<'>{p, z, z') = A<'> zlj^^.ipz) + B<'> zlY^^^Jpz), 



(F.24) 



where J„ and Yn are Bessel functions of the first and second kinds, A^'^ and i?^'^ are coefficients 
to be determined by boundary conditions, and p is the analog of Xp defined hj p = y^p^p^^. Note 
that this differs from (F. 9) since there is no explicit bulk mass dependence. In ( F.24[ ) the bulk 
masses enter only in the order of the Bessel functions as (c ± |). 



The matching boundary conditions at z = z' are given by (F. 10) and (F. 11) modified by a 



factor of {R/z')^ from (F.23), 



d5F>(z')-d,F<(z' 



5^± 



F>{z')-F<iz') 



t{R/z')-\ 
0. 



(F.25) 
(F.26) 



The chiral boundary conditions are the same as in the flat case, (F. 12) and (F. 13) with the 
appropriate insertion of (F.21). 

We may now solve for the A and B coefficients. It is useful to write these in terms of common 
factors that appear in their expressions. To this end, let us define the prefactors 



ITT 



1 



«L 



2R^ S-{pR,pR') 



ITT 



2R^ S+{pR,pR') 



(F.27) 
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A^< = 




_._{pR)St{pz\pIi) 


Al< = 


-anz'K^ ipR)S+{pz',pR') 


B^< = 




_^{pR)~St{pz\pR!) 


Bf< = 


aRz'lj^^{pR)S+ipz',pR') 


A^< = 




^._{pR)S-{pz',pE!) 


A^< = 


~aRz'K+\ {pR)S-{pz',pR') 


B^< = 


aiz 2 


_r{pR)S-{pz\pE!) 


B^< = 


aRz'iJ^^i (pR) S~{pz\pR!) 



Table 2: Left-handed RS fermion propagator coefficients: the z > z' coefiicients are obtained by 
swapping R R' va. the arguments of the functions, leaving the aL,R constant. 



and a set of antisymmetric functions 



c±i(y)^c±i(a;) 



(F.28) 
(F.29) 



With these definitions the coefficients for the left- and right-handed F functions are given in 
Table |2l The Ff'^ functions may thus be written out succinctly for 2; < z' as 



F^< = aL {zz'f^ S+ {pz',pR') S-{pR,pz) (F.30) 

Ff< = ai {zz'fl'^ S- {pz',pR') S;{pR,pz) (F.31) 

F^< = an {zz'f' S+ {pz',pR') SUpR,Pz) (F.32) 

F^< = aR {zz' f^ S; {pz'.pR!) StipR^pz) (F.33) 

The expressions for z > z' are obtained by making the replacement {R ^ R'} in the arguments 

of the Sc functions. We now use the notation in (F. 7) and drop the <, > superscripts. From these 
the fermion Green's function can be obtained trivially from the analog of (F. 6), 

A{p,z,z') = VF{p,z,z') -- 



a^Pf^F^ D^F+ 



D, 



± 9; 



+ 



(F.34) 



Note that in the UV limit {xp 1/-R) the Bessel functions reduce to phase-shifted trigonometric 
functions so that we indeed recover the fiat 5D propagators. 



F. 3 Euclidean warped 5D fermion propagator 

Finally, it is convenient to write the Wick-rotated form of the fermion propagators since these 
will provide the relevant Feynman rules in loop diagrams such as /i — )■ 67. We shall write out the 
scalar F functions in a convenient form that we use throughout the rest of this document. The 
derivation is identical to that outhned above with the replacement = — p|; (i. e. d = iDe) in the 
Green's function equation so that we shall simply state the results. The Euclidean scalar functions 
are written in terms of the modified Bessel functions / and K which behave like exponentials in 
the UV. Let us define the auxiliary functions 

Sc{x±,x'j^) = Ic±l/2{x)Kc±l/2{x') - Ic±l/2{x')Kc±l/2{x) (F. 35) 

Scix±,x'^) = Ic±l/2{x)Kc^i/2{x') - IcTl/2{x')Kc±l/2{x) (F. 36) 

Tc(x±,x'^) = Ic±l/2ix)Kc^i/2{x') + IcTl/2ix')Kc±l/2{x). (F. 37) 
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Since we would like to write dimensionless loop integrals, let us define the dimensionless variables 
y = UbB! and x = ksz, which are the natural quantities which appear as arguments of the Bessel 
functions. We write the warp factor as w = {R/ R'). It is convenient to puli out overall factors to 
write the F functions as 



F±{kE,z,z')=tw-^R'F^^y. 
The Euclidean scalar functions for x > x' (i. e. z > z') are given by 

ixx')'/' S,,{x.,y^)S,,{x'_,wy.) _ {xx')'/' n,ix+,y^)n,{x'^,wy^) 



F!l 



y' 



(XX 



')^/^T,,{x_,y+)n,{x'_,wy^) 



(XX 



(F.38) 

(F.39) 
(F.40) 



Scn{y+,wy+) + 7/5 Scj,iy+,wy+) 

The functions for x < x' are given by replacing x -H- x' in the above formulas. With these 
definitions the Euclidean fermion propagator given by the analog of (F.34), 



^ ' ' ^ y \a^y^F_ yD_F+) 



D+ = ± dr- 



+ 



X 



(F.41) 



G Finiteness of the brane-localized neutral Higgs diagram 



As explained in Section |6.4[ the finiteness of the one-loop result and logarithmic divergence at 
two-loop order becomes opaque to nai've 5D power counting arguments when the Higgs is brane- 
localized. Additional cancellations of leading-order terms in loop momentum are required to 
sensibly interpolate between the superficial degree of divergence of the bulk and brane-localized 
scenarios. For the charged Higgs this cancellation mechanism came from an insertion, which 
led to an additional l/k"^ factor relative to the bulk field. Here we shall elucidate the finiteness of 
the single-mass-insertion brane-localized neutral scalar loop. 

At one-loop order this finiteness can be seen explicitly by the cancellation between the neutral 
Higgs and the neutral Goldstone. However, there is an additional chiral cancellation that occurs 
between the two diagrams associated a single intermediate neutral boson. Indeed, because the 
Higgs and neutral Goldstone do not appear to completely cancel at two-loop order, this additional 



cancellation is necessary for the power-counting arguments given in Section |6.7 

We highlight this cancellation in two ways. The pure momentum space calculation highlights 
the role of the chiral boundary conditions, while the mixed position/momentum space calculation 
shows an explicit cancellation while including the fuU scalar structure the amplitude. 



G.l Momentum space 

Here we shall see that 4D Lorentz invariance combined with the chiral boundary conditions forces 
the UV divergence of the two diagrams in Fig. |8]to cancel. 

We first note that the propagators to the photon vertex each have an endpoint in the bulk. This 
implies that the leading-order contributions to these propagators in the UV limit are proportional 
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to the uncompactified flat-space 5D propagators, 



^A,, A^J ^ - kl [k.a'^ -th) k^-kl ' ^"^-'^ 

where we have written A^^ to mean the propagation of a left-handed Weyl spinor x i^to a right- 
handed spinor ip. The terms along the diagonal come from /^57^ and represent the chirahty-flipping 
part of the propagator. The boundary conditions require the wrong-chirahty modes, the SU(2) 
doublet iIjl and SU(2) singlet xRi to vanish on the IR brane. Thus, the fermion may propagate 
to the wrong-chirality spinor in the bulk only if it propagates back to the correct-chirality spinor 
when it returns to the brane. For an internal left-handed Weyl fermion xl-, the portion of the 
amplitude coming from the photon emission takes the form 

A^^a^A^^ + A^^a^A^^ ~ (A;„a") [k^a^) + (^5) (G.2) 

Combining with the analogous expression for a right-handed Weyl fermion in the loop, the relevant 
part of the photon emission amplitude can be written as 

(P - klY ' ^ ^ ^ 

where these terms correspond to a fermion of the correct and incorrect chirality propagating into 
the brane. The second term can be simplified using 

{k^Y f ~k'^ 

which can be confirmed by Wick rotating both sides, k"^ — )• —k"^, and performing the dfcs integral 
explicitly. Now it is easy to see that the divergent contributions from the diagrams in Fig. [8] 
cancel. The boundary conditions force brane-to-brane propagators to go like ^ with no 7^ part. 
Thus we may write the internal fermion structure of the amplitudes as 

Mia) + Mib) ~ ^ ih'"^ - ^V) + ih'"^ - ^'7'') 1^ = 0. (G.5) 

The key minus sign between the two terms in the photon emission comes from the chiral boundary 
conditions that force the second term to pick up the relative sign between the two diagonal blocks 
of 7^. 



Let us remark that it is crucial that the denominator in (G.4) contains exactly two propagators 
or else the equality would not hold. One might be concerned that the brane-to-brane propagator 
should also contribute an additional factor of (fc^ — /c|) to the denominator (the k^'-f^ term vanishes 
in the numerator from boundary conditions). Such a factor is indeed present in the fuU calculation, 
but because 5D Lorentz invariance is broken on the brane, k^ is not conserved there and this factor 
actually includes a different, uncorrelated fifth momentum component, ^5, which can be taken the 
be independent of the dk^ integral. This is a manifestation of the principles in Appendix [A| As a 
check, one can perform the dk^ integral for this brane-to-brane propagator and obtain the same 



\k\ UV behavior found in the careful derivation performed in Appendix F.l 
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G. 2 Position/momentum space 



In Appendix |F.1| we derived the flat-space bulk fermion propagator, 

A{p,x„4) = {p-^^% + m) - '^^ - ^^^'^ + ^'^("^^ ^"^^^ ' + ^'^^^ , (G.6) 

2xp sin XpL 

where the zero mode chirality is given by X = {L, R} with p{L) = +1 and p{R) = —1. We then 
argued at the end of Appendix F.2 that the propagators in a warped extra dimension reduce to 
this case up to overall phases. Thus we expect the amphtudes to have the same UV behavior up 
to finite factors. The relevant flat-space one-loop diagrams contributing to the operator (4.1) are 
shown in Fig. [Sj We start with Fig. 8a and assume that the decay is from /iL to cr. The loop 
propagators with {x5,x'^) = {L, z), {z, L) and {L, L) can be written as 

A/,/ r ^ .}^' cos Xk' z -i-f^Xk' sin Xk' z , . 

A{k,L,z) = -i : Pr (G.7) 

Xk' sinXk'L 

T\ . f(: cos XkZ + i'y^XkSmXkZ , . 

A{k,z,L) = -i : Pr (G.8) 

Xk sinXkL 

A{k,L,L) = -i : -Pr, (G.9) 

Xk sinXkL 

where k' = k + q. We have used the chiral boundary conditions to simphfy A (/c, L, L). Since we 
are interested in the UV behavior we have dropped the terms proportional to the bulk mass m 
from the internal propagators because these are finite. Combining the propagators together and 
doing the same calculation for Fig. |8bl the amplitudes become 



\Ai^ f ^ -( Y^f{k,z)+XkXk'l^g{k,z)\]/^cotXkL 

= / jT^dz u{p)\ ^ \ u{p) G.IO 

^''^ J (27r)4 1^ XkXk' [{p + ky - mjj\ J 



Xk 



d^k ^ _^^,J' cotXk'L \H'l^^f{k,z) + XkXk'l''9{k,z 
where we have written 



^ ' I (2vr)4 Xk' [ XkXk' [{p + ky - mjj] 



/(fc,.) = -^^^^^^±^^^l^^^ (G.12) 
sm Xk+gL sm XkL 

smXk+qLsmXkL 

Note that all of the z dependence is manifestly contained in sines and cosines. Further we have 
neglected the flavor-dependence of the Xk factors since these also come from the bulk masses via 
(F. 9) and are negligible in the U V. 

Upon Wick rotation the trigonometric functions become hyperbolic functions which are expo- 
nentials in the Euclidean momentum, 

cos XkZ cosh{xkEz) = ^ {e^^E^ + e"'^*^^) (G. 14) 

sinxfc-2 ismh{xkEz) = - (e^'=s^ - e'^'^B^) . (G. 15) 
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We may now replace the trigonometric functions with the appropriate Euclidean exponentials. 
Since we are concerned with the UV behavior, we may drop terms which are exponentially sup- 
pressed for large k over the entire range of z. The remaining terms are simple exponentials and 



can be integrated over the interval. One finds that the trigonometric terms in (G. 10) and (G. 11) 
yield the expression 



XkE+q ~l~ XkE Xk+q Xk 



(G.16) 



where on the right we have reversed our Wick rotation to obtain a Minkowski space expression for 
the terms which are not exponentially suppressed in Euchdean momentum. After doing this, the 



leading order term in cot xL in (G. 10) and (G. 11) equals i ^ and the terms in the braces become 



- Xk+qXk Y 



. XkXk+q {Xk + Xk+q) [{p + kf - 



m 



H\ 



(G.17) 



which gives the numerator of (G.5). 

In terms of these quantities the potentially divergent amphtudes can be written as 



(a) 



1 



-u{p) 



+ 



d^k 

(2vr)4 {xk+q + Xk)[ij> + ^? - m]j\ ' [ Xk+q 

d k 1 - { \ f ^ - /- 

{Xk+q + Xk)[{p + kr-mlCP' V Vk~ 1^ 



+ 



Xk ) 
^-Y\u{p + q), 



(G.18) 
(G.19) 



therefore these two terms cancel each other in the UV and the operator (4.1) is finite. 



Higher mass insertions do not spoil this cancellation since these are associated with internal 
brane-to-brane propagators whose UV limit goes like A(A;) ~ t/Xk- The chiral structure of 



the effective operator (4.1) requires that only diagrams with an odd number of mass insertions 
contribute. Using the UV limit A(A;)^ —t- 1 one notes that the divergence structure reduces to the 
case above. 
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